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Institute for Nuclear and Energy Technologies 


Structure and Laboratories of the Institute of Nuclear and 


Energy Technologies 


Walter Tromm 


The Institute for Nuclear and Energy Technol- 
ogies (IKET) is situated with its offices and re- 
search laboratories in Campus North of Karls- 
ruhe Institute of Technology (KIT). 


The research of IKET focuses on conversion 
from thermal power to electric power for future 
power plants without CO2 emission, on hydro- 
gen technologies for energy storage and on 
safety. Traditionally, nuclear power plants, for 
which IKET concentrates on their safety fea- 
tures and on methods to mitigate severe acci- 
dents, were in the focus of investigation. IKET 
uses its competence, gained from a profound 
experience in numerical simulation and in de- 
sign and operation of large-scale experimental 
facilities to apply this knowledge more and 
more in the area of renewable energies. 


The combination of science and technology 
with education and training is a systematic ap- 
proach at KIT, and IKET is contributing accord- 
ingly to courses in mechanical engineering, su- 
pervises several bachelor and master theses 
each year and coordinates master programs in 
energy technologies. Compact courses on en- 
ergy technologies are given also in executive 
master programs and in the Framatome Pro- 
fessional School, which is funded by industry 
and managed by IKET. 


Resources 


The employees of IKET contribute in the re- 
search field energy of the Helmholtz Associa- 
tion (HGF) to the research-programs Nuclear 
Waste Management, Safety and Radiation Re- 
search (NUSAFE), Fusion (FUSION), Renew- 
able Energies (EE) as well as Storage and 
Cross-linked Infrastructures (SCI), with a focus 


on NUSAFE. By the end of 2019, around 75 
Scientists, engineers and technicians have 
been working at IKET on this wide range of 
CO»-free technologies for energy conversion. 
Fig. 1 illustrates that more than half of the IKET 
employees were still working for the NUSAFE 
program, despite the German phase-out deci- 
sion, but in accordance with the German en- 
ergy politics and the 7!" Energy Research Pro- 
gram of the German Federal Government. 
Around two third of the employees were 
funded in 2019 by the Helmholtz Association 
(HGF), the others by third party funds of the 
European Commission, by industry, by Ger- 
man ministries or by other research funds. 
Doctoral students as well as students of the 
Baden-Wuerttemberg Cooperative State Uni- 
versity (DHBW) were filling around 2096 of 
these positions at IKET. In addition, students 
perform their bachelor or master theses or 
spend an internship in the research laborato- 
ries of IKET. 


2019 


BFUSION @NUSAFE BEE nSCI 


Figure 1: Assignment of IKET personnel to HGF pro- 
grams. 
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An overview of the structure of IKET is given 
by the organization chart, Fig. 2. Each working 
group is acting in close collaboration with the 
other groups, particularly in its own research 
program. 
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Figure 2: Organization chart of the Institute for Nuclear 
and Energy Technologies 


Working groups on Severe Accident Research, 
and on Transmutation as well as the Frama- 
tome Professional School have been concen- 
trating in 2019 primarily on nuclear applica- 
tions, whereas the Karlsruhe Liquid Metal 
Laboratory (KALLA) and the Hydrogen group 
were addressing nuclear as well as renewable 
energy technologies. The working group on 
Magneto-Hydrodynamics is primarily working 
on nuclear fusion applications, whereas the 
working group on Energy and Process Engi- 
neering is rather concentrating on geothermal 
energies. The working group Multiphase Flows 
is mainly experimentally investigating techni- 
cally relevant, complex flow phenomena. The 
field of activity currently comprises mainly the 
complex structure of the COSMOS-H high- 
pressure loop. The working group on Accident 
Management Systems continued to extend the 
application of its simulation models to several 
critical infrastructures. Thus, the institute co- 
vers a wide field of different energy technolo- 
gies, and the share of its personnel resources 
on the different research topics reflects the re- 
quirements of the Helmholtz Association's en- 
ergy research field. 


They are all supported by a joint infrastructure, 
comprising a metal workshop, manufacturing 
urgent test components, a welding shop, and 
an electromechanical workshop. Other tasks of 
the infrastructure include the IT-administration, 
business administration and public website of 
the IKET. The Infrastructure team is active as 
well in education and training activities. 


Every year, at least four students of the Baden- 
Württemberg Cooperative State University are 
employed by IKET, managed by the Infrastruc- 
ture group, to work with the research teams as 
part of their educational program. 


Working in a research area with industrial ap- 
plication IKET practices a Quality Management 
System appropriate to EN ISO 9001. 
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International Research Projects 


Most studies performed at IKET are embedded 
in international research projects. A few ex- 
traordinary highlights for 2019 are given below: 


Highlight: Hydrogen group: 
HyTunnel-CS 


The Hydrogen group is strongly involved in the 
HyTunnel-CS project. The aim of the HyTun- 
nel-CS project is to perform pre-normative re- 
search for safety of hydrogen driven vehicles 
and transport through tunnels and similar con- 
fined spaces. The main ambition is to facilitate 
hydrogen vehicles entering underground traffic 
systems at risk below or the same as for fossil 
fuel transport. The specific objectives are: crit- 
ical analysis of effectiveness of conventional 
safety measures for hydrogen incidents; gen- 
eration of unique experimental data using the 
best European hydrogen safety research facil- 
ities and three real tunnels; understanding of 
relevant physics to underpin the advancement 
of hydrogen safety engineering; innovative ex- 
plosion and fire prevention and mitigation strat- 
egies; new validated CFD and FE models for 
consequences analysis; new engineering cor- 
relations for novel quantitative risk assessment 
methodology tailored for tunnels and under- 
ground parking; harmonised recommenda- 
tions for intervention strategies and tactics for 
first responders; recommendations for inher- 
ently safer use of hydrogen vehicles in under- 
ground transportation systems; recommenda- 
tions for RCS. 


The objectives will be achieved by conducting 
inter-disciplinary and inter-sectoral research by 
a carefully built consortium of academia, emer- 
gency services, research and standard devel- 
opment organisations, who have extensive ex- 
perience from work on hydrogen safety and 
safety in tunnels and other confined spaces. 
The complementarities and synergies of theo- 
retical, numerical and experimental research 
will be used to close knowledge gaps and re- 
solve technological bottlenecks in safe use of 


hydrogen in confined spaces. The project out- 
comes will be reflected in appropriate recom- 
mendations, models and correlations. HyTun- 
nel-CS will reduce over-conservatism, 
increase efficiency of installed safety equip- 
ment and systems to save costs of under- 
ground traffic systems. 


The 1°! progress meeting was hosted in Karls- 
ruhe by the Hydrogen working group on the 
11-131 September 2019. 


Figure 3: HyTunnel-CS consortium and SAB members at 
the 1st progress meeting in Karlsruhe, 13th September 
2019. 


Highlight UNF: CONFIDENCE 


First of January 2017, the CONFIDENCE 
(COping with uNcertainties For Improved mod- 
elling and DEcision making in Nuclear emer- 
genCiEs) project started as Task 9.1 of the 
H2020 project CONCERT. The project focuses 
on identifying and reducing uncertainties in the 
release and post-release phases of an emer- 
gency. The latter includes the transition be- 
tween the short-term post-release and recov- 
ery phases (e.g. the first year(s)). The group 
emergency management (UNF) coordinated 
CONFIDENCE. The scientific challenges were 
addressed though six interlinked work pack- 
ages (WP), with an additional E&T WP. UNF 
also contributed to several work packages and 
coordinated WP6. 


Besides the improved Multi Criteria Decision 
Analysis tool described in this report, UNF was 
heavily engaged in enhancing the decision 
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support system JRODOS for nuclear and radi- 
ological emergencies to deal with uncertain in- 
put of weather and source term. In this respect, 
JRODOS was extended by an ensemble gen- 
erator and an ensemble visualisation engine 
(ensemble approach means realisation of a 
state by many simulations with slightly different 
initial conditions). A typical application case 
consists of an ensemble of 20 weather files 
and 3 source terms (realistic, worst and opti- 
mistic case) with their probabilities of occur- 
rence. The resulting 60 scenarios are pro- 
cessed in one project and stored individually. 
In a second step and evaluation module allows 
visualising probability functions of interesting 
results such as areas of early countermeas- 
ures. The percentiles represent probability that 
a grid point is affected by a particular counter- 
measure. The 50% percentile means that in 
our example with 60 cases, 30 result in a coun- 
termeasure at that location. Visualisation is still 
an issue that needs further discussion with de- 
cision makers. Colour tables of the visualisa- 
tion can be adapted by the user 


Figure 4: Map from the JRODOS ensemble processor for 
ensemble runs with blue colour table and added contour 

plots separating areas with fraction borders of 1, 0.8, 0.5 

and 0.2 


Future research will concentrate on the testing 
of the ensemble approach in particular related 
to operational applications, communicating to 
decision makers and lay people as well as 
comparing different weather providers with 
each other. There is one important question 


not resolved: Do the numerical weather predic- 
tion ensemble really represent the uncertainty 
of the weather? 


Highlight KALLA working group: 


Liquid metals are tested as heat transfer fluids 
in concentrating solar power systems in the 
SOMMER test facility (Fig. 1) at Karlsruhe Lig- 
uid Metal Laboratory KALLA. The excellent 
heat transfer properties of liquid metals allow 
building more compact and efficient receivers. 
The thermal receiver tubes in the focal point, 
which are heated by the concentrated solar 
power, can be cooled more efficiently than with 
the currently used molten salt mixture. In the 
SOMMER loop, the liquid lead-bismuth eutec- 
tic (LBE) is pumped through the spiral thermal 
receiver and is afterwards cooled down in an 
air cooler to enter the pump again. In 2019, ex- 
periments under sunlight conditions with heat 
flux densities of up to 4 MW/m? have been suc- 
cessfully performed in the SOMMER test facil- 
ity. It was demonstrated that liquid metal can 
be handled under varying conditions and can 
efficiently cool the receiver pipes at extremely 
high heat flux densities. 


Figure 5: SOMMER test facility; left: parabolic mirror; 
right: liquid metal loop with thermal receiver in the focal 
point (Foto: F. Fellmoser) 


Institute for Nuclear and Energy Technologies 


Hydrogen from Natural gas without 
CO2 Emissions 


c^ 
\ 


iret 


Figure 6: Methane pyrolysis by means of a bubble col- 
umn reactor allows for the climate-friendly use of fossil 
natural gas. 


In the current energy debate, hydrogen is in- 
creasingly considered a key to the success of 
the energy transition. Hydrogen is a clean 
source of energy and can be used for the pro- 
duction of electricity and heat and in the mobil- 
ity sector or in industrial processes, such as for 
the production of steel. Experts of the Interna- 
tional Energy Agency IEA calculated that ad- 
mixture of 2096 of hydrogen in the European 
gas grid could reduce CO2 emissions by 60 
million tons per year. This corresponds to 
about the amount emitted by Denmark in a 
whole year. 


Direct thermal cracking of methane and other 
hydrocarbons can be a way to produce hydro- 
gen from natural gas without direct CO2 emis- 
sions. For this purpose, the KALLA team, in co- 
operation with the Institute for Advanced 
Sustainability Studies e.V. in Potsdam, started 
developing a process for the continuous de- 
composition of methane in a bubble column re- 
actor filled with liquid metal to produce hydro- 
gen and solid carbon some years ago. They 
continued research on the technology over the 
last years and have entered into a joint project 
with the industry partner WintershallDea in 
2019, to further advance this technology to- 
wards use in industry. 


Highlight SAR-group: JIMEC 
experiments for the EU-Project ESFR- 
SMART 


To increase public acceptance of nuclear 
power in Europe and secure its future role, the 
significantly higher safety of new reactors com- 
pared to traditional reactors has to be demon- 
strated. The 4-year ESFR-SMART project (Eu- 
ropean Sodium Fast Reactor Safety Measures 
Assessment and Research Tools), launched in 
September 2017 aims to enhance the safety of 
Generation-IV Sodium Fast Reactors (SFR), in 
particular the commercial-size European So- 
dium Fast Reactor (ESFR). 


In course of a less probable accident that a co- 
rium pool builds up in the core region, dedi- 
cated corium relocation measures with Corium 
Transfer Tubes will minimise the core melt 
mass in the core region to avoid recriticality. 


A great safety concern in this case is the im- 
pingement of a high-temperature and high-mo- 
mentum corium jet on the metallic core catcher 
material. Although both oxide corium and me- 
tallic corium can be relocated, metallic corium 
is considered to be more crucial on the core 
catcher ablation process due to its low solidify- 
ing temperature and therefore the less poten- 
tial to build up a protecting crust layer. 


The JIMEC (Jet Impingment on Metallic Core- 
catcher) tests performed by the Severe Acci- 
dent Group (SAR) of IKET provided prototypi- 
cal experimental data to quantify the risk. 
JIMEC tests are designed firstly to create a 
metallic jet with stainless steel composition 
and with designed superheat and fall height; 
and secondly to record the ablation process of 
a 40 cm thick substrate of the same material 
composition as the jet. The mass of the metal- 
lic melt is 1000 kg, and the jet velocity at the 
substrate surface is about 5 m/s. 


The prototypical test data provided insights to 
a wide application range within and also be- 
yond the nuclear safety research. 


The main conclusion of the experiments is that 
the ablation process slows down after the melt 
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jet mass was firstly retained in the molten pool 
inside the ablated pit. The ablation velocity be- 
fore and after this event was constant with an 
ablation velocity of less than 2 cm/s before and 
slowed down afterwards to less than 1 cm/s. 
This effect could be called "pool effect", which 
means that the deep melt pool prevents the di- 
rect energetic contact of the melt at the pit bot- 
tom. 


0 sec: outburst of metallic jet 


22 sec: pool effect in the sub- 72 sec: direction of oxide jet after- 
strate wards 


im u on 


Figure 7: shows the impingement of a 40 mm diameter steel jet on a steel substrate: on the top, the jet has direct 
contact with the substrate during the initial period, while at the bottom, the jet falls into a molten pool and reduces its 
ablation effect after a certain depth reached in the ablated pit 
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Magnetohydrodynamics for liquid-metal blankets 


L. Bühler, H.-J. Brinkmann, C. Courtessole, V. Klüber, C. Köhly, C. Mistrangelo 


Introduction 


Liquid metal blanket concepts foreseen in a 
DEMOnstration nuclear fusion reactor or as 
test blanket modules (TBM) in the experimen- 
tal reactor ITER use the lead lithium alloy PbLi 
as breeder, neutron multiplier, neutron shield 
and as coolant to remove the volumetrically 
generated heat. The movement of the electri- 
cally conducting fluid in the strong magnetic 
field induces electric currents responsible for 
strong Lorentz forces, high magnetohydrody- 
namic (MHD) pressure drop with implications 
on heat and mass transfer. The physical pa- 
rameters characterizing the MHD flow are the 
Hartmann number and the Reynolds number, 


Ha = BL |S Rez 
pv v 


where B, L and uo denote magnetic field, typi- 
cal length scale and velocity. The fluid proper- 
ties are the density p, kinematic viscosity v and 
electric conductivity o. Ha? und Re stand for 
the ratio of electromagnetic to viscous forces 
and inertia to viscous forces, respectively. For 
problems, where buoyancy forces are im- 
portant, the latter ones are quantified by the 
Grashof number 


ATI? 
NEU : 
v 


’ 


where AT is a typical temperature difference, g 
stands for gravity and £ is the thermal expan- 
sion coefficient of the liquid metal. 


The MHD work at IKET KIT supports blanket 
engineering activities of the EUROfusion con- 
sortium and the test blanket program of Fusion 
For Energy (F4E) for ITER. 


Numerical simulations for liquid metal 
MHD flows 


Optimum grids for MHD simulations 


Phenomena imposed by high transverse mag- 
netic fields result in strong electric current den- 
sity and steep gradients of velocity across ex- 
tremely thin boundary layers. The later have to 
be spatially well resolved in numerical simula- 
tions since they determine essentially the flow 
in the entire duct. In contrast, the flow in the 
core shows only moderate or small gradients 
so that this part of the flow could be simulated 
with coarser grids. Computational grids based 
on cubic cells are favorable for finite-volume 
calculations of MHD flows in rectangular ducts. 
As shown in Figure 1, by stretching and com- 
pressing the grid independently in both trans- 
versal directions, block structured cells adapt 
perfectly to the duct geometry, while preserv- 
ing a structural order, that benefits to numerical 
consistency. 


Figure 1: Computational grid for cross section of a rec- 
tangular duct MHD flow. 


For more complex geometries block-structured 
grids might not be an optimal choice. Ongoing 
investigations therefore address the applica- 
tion of unstructured grids for simulations with 
an in-house-developed MHD code based on 
the finite-volume library OpenFOAM [1]. One 
example discussed in the following is the well- 
known MHD flow in an electrically insulating 
circular pipe. 
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rectangular 


polyhedral block-structure 


Figure 2: Three grid types for pipe flow and their re- 
spective impact on axial MHD velocity profiles (lower 
row). 


Figure 2 shows different types of grids together 
with their resulting velocity profiles. All grids 
use hexahedral prisms in order to resolve thin 
boundary layers, while using rectangular, poly- 
hedral or a block structure in the core. Rectan- 
gular and polyhedral based meshes show cor- 
rect MHD velocity profiles for circular pipe flow. 
We observe only minor ripples at positions with 
locally unstructured cells. The block-structured 
configuration (standard O-grid), however, ex- 
hibits large singularities emerging from the cor- 
ners of the quadratic central core grid. These 
strong perturbations propagate along grid lines 
into the fluid. This seems to be an intrinsic 
problem in finite volume calculations as re- 
ported in [1] and [2]. According to [3], locations 
of slope discontinuity, which are characteristic 
for the block-structured O-grids, lead to zero 
order errors, independent of the mesh refine- 
ment level. Computations of MHD flows and in 
particular, the strong coupling of hydrodynamic 
and electrodynamic phenomena amplifies this 
effect additionally. Consequently, rectangular 
or polyhedral grids should be the preferred 
choice for simulations of MHD flows in complex 
geometries. 


3D MHD flows in fusion blankets 


The dual coolant lead lithium (DCLL) blanket is 
one concept investigated in the frame of EU- 
ROfusion [4]. In the proposed design, PbLi 
serves as breeder and as coolant. Since the 
liquid metal flows at sufficiently large velocity 


for efficient heat removal, pressure losses re- 
lated to the occurrence of electromagnetic 
forces become significant. Magnetohydro-dy- 
namic effects are coupled with thermal phe- 
nomena caused by the occurrence of temper- 
ature gradients in the fluid due to the presence 
of a volumetric neutronic heat load. 


In order to reduce MHD pressure drop it has 
been proposed to use insulating flow channel 
inserts (FCl) to electrically decouple the PbLi 
flow from the well-conducting walls. Sandwich- 
type FCls are considered as a feasible tech- 
nical option [4] [5] [6] and their electrical con- 
ductivity (of the inner protecting sheet) has 
been taken into account in numerical simula- 
tions. 


N 
~ 

Q 

a 


tory rad, x 


A 


Figure 3: (a) Design of DCLL blanket module. (b) Model 
geometry for the numerical study. 


Sa 


From Figure 3, it can be noticed that the main 
characteristic geometric elements of a module 
area vertical rectangular ducts. Initial numeri- 
cal simulations for pressure driven flow have 
been complemented by new 3D calculations 
for mixed convection of PbLi in the presence of 
a strong toroidal magnetic field. The flow is 
driven by an external pressure gradient and by 
buoyancy associated with the volumetric heat- 
ing caused by neutrons. A model problem for 
mixed magneto-convective flow is considered, 
where a uniform volumetric heat source q in 
the fluid is assumed. The characteristic tem- 
perature difference is related to this heat 
source as AT = qI?/k, where k stands for the 
thermal conductivity of the liquid metal. At the 
entrance, the MHD flow is fully developed and 
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isothermal. In the inlet duct and in the entire 
blanket module, the fluid heats up along the 
flow path as a result of the volumetric heating. 


Figure 4 shows results for mixed convection 
flow at Re= 8596, Ha = 500 and Gr = 10°. The 
temperature distribution on the middle plane at 
y = 0 shows that a high-temperature region is 
present on the top of the module and in the up- 
per part of the “downward” duct. Streamlines 
and contours of velocity magnitude depicted in 
the figure indicate the presence in the rear 
channel of a descending jet that moves down- 
stream along the internal vertical wall that sep- 
arates the two poloidal channels. Here the fluid 
is colder and therefore, due to its higher den- 
sity, it tends to flow downwards. The downward 
jet is unstable and this may affect mixing and 
the heat transfer. 
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Figure 4: Results for Re= 8596 Ha = 500 and Gr = 
10? on the plane y = 0; (a) contours of temperature and 
velocity streamlines; (b) contours of velocity magnitude. 


Velocity profiles for forced and mixed convec- 
tion at Ha = 500 and Re = 8596 are shown in 
Figure 5 for Gr = 0 and Gr = 10?, respectively. 
The intensity of the downward jet flow at the 
internal wall can be seen. More details can be 
found in [7] [8]. Results for magneto convection 
in geometries related to water cooled lead lith- 
ium WCLL blankets are described in [9]. 
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Figure 5: Velocity distributions at different positions 
along the flow path for MHD flow at Ha = 500, Re = 
8596 for (a) isothermal flow Gr = 0, (b) magneto con- 
vection at Gr = 10°. 


Experimental study of liquid metal 
MHD flows at the entrance of FCls 


MHD pressure drop in ducts of fusion blankets 
results from the interaction of flow-induced 
electric currents with the plasma-confining 
magnetic field. The highest pressure drop in 
the flow is observed, when currents shortcut 
along the thick conducting blanket walls. For 
reduction of pressure drop, electrically insulat- 
ing flow channel inserts (FCls) are foreseen for 
electrically decoupling the liquid metal flow 
from the well-conducting channel walls. The 
preferred option for near-term applications in 
fusion blankets are sandwich-type FCls, where 
the insulating ceramic material is protected 
from direct contact with the liquid metal by thin 
sheets of steel. For experimental investiga- 
tions of 3D effects at the entrance into FCls 
and at junctions between two FCls, a test sec- 
tion has been manufactured and experiments 
have been performed in the MEKKA facility 
[10] using NaK as model fluid. The fabrication 
of FCls for circular pipes required development 
of new fabrication strategies [6], since the tech- 
nology proposed for rectangular ducts [5] did 
not apply. The experimental study shows the 
benefits of FCls for pressure drop reduction in 
fully developed flows as predicted by theoreti- 
cal analyses [11] and confirmed by experi- 
ments [12]. 
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Details of the experimental test section are 
shown in Figure 6. The thickness tpc = 
0.5mm of the protection sheets has been sug- 
gested e.g. in [13]. The coordinate x = 0 cor- 
responds to the position of the pressure tap 
that is located immediately in front of the FCI. 
As shown in Figure 6 and Figure 7, the FCI 
starts at x = 5mm and the insulation inside the 
FCI at x; 2 10mm. 


3D perturbations that occur when the flow en- 
ters the magnet decay quickly along the flow 
path as shown in previous investigations [14]. 
Therefore, the MHD pipe flow is fully devel- 
oped before the fluid approaches the FCI. In 
the experiments, pressure and electric poten- 
tial differences have been recorded and some 
typical results are discussed in the following. 
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Figure 6: Geometry and principle sketch of current 
paths at the entrance of a FCI. The insulation inside the 
FCI starts at x; = 10mm. 
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Figure 7: Sketch of FCI position in pipe. The insulating 
layer inside the FCI (not visible in the sketch) begins at 
a downstream position of x; = 10mm. 
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In experiments it is possible to observe a po- 
tential difference between both sides of the 
pipe Ad= (a = -2) - $ (a = 2). This 
quantity is used in the following to investigate 
the influence of 3D MHD effects at the en- 
trance of the FCI. When potential results are 
scaled by characteristic values uoBL, all curves 
collapse onto a single line as displayed in Fig- 
ure 8, confirming the universal scaling law for 
potential over a wide range of Hartmann num- 
bers investigated. Potential values are propor- 
tional to the strength of the applied magnetic 
field and the fluid velocity. 


Far upstream, A@(x/L < —2) shows fully es- 
tablished conditions with uniform values along 
x in good agreement with theoretical predic- 
tions for fully developed MHD pipe flow [15]. 
Approaching the FCI, the potential starts mon- 
otonically decreasing. At downstream posi- 
tions, where current paths into the wall are in- 
terrupted by the insulation of the FCI, values of 
AQ (x/L) decay to zero. We observe that 3D ef- 
fects, caused by recirculating currents inside 
the wall near the entrance of the FCI, extend 
upstream over more than 2L and downstream 
up to 5L. The currents flowing in the down- 
stream part of the wall are separated from 
those inside the fluid by the insulation in the 
FCI. The fact that Ab vanishes as x/L > 5 con- 
firms full functionality of the ceramic insulation 
inside the sandwich FCI. 
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Figure 8: Variation of nondimensional transverse poten- 
tial difference near the entrance of the FCI. 
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Pressure differences along the pipe and along 
the FCI entrance region have been measured 
for various values of flow rate and for different 
strengths of the magnetic field. The reference 
value of pressure has been placed at position 
x = xj, i.e. p(x;) = 0. When scaled with char- 
acteristic quantities ou,B?L, all results come 
close to a single line as shown in Figure 9. We 
Observe a fast decay of pressure in the bare 
pipe without insulation and a slower decrease 
of pressure in the insulated region. 
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Figure 9: Nondimensional pressure along the axis for 
Ha = 2000 and different. 


The pressure gradients in the bare pipe and in 
the part with FCI agree well with theoretical 
predictions [16] as indicated by the black solid 
lines in Figure 9. The FCI reduces the pressure 
gradient by a factor of about 13.3 compared to 
the value in the non-insulated pipe. Neverthe- 
less, recirculating additional currents near the 
entrance of the FCI give rise to extra Lorentz 
forces that create some additional pressure 
drop Ap3p caused by 3D effects. Even if the ab- 
solute value of App appears acceptable one 
should keep in mind that it corresponds to a 
pressure drop in an ideal long FCI over more 
than 4 characteristic lengths L. Since similar 
pressure losses occur also at the exit of the 
FCI or at gaps between two FCls [12], these 
3D effects reduce the efficiency of FCls by an 
amount that is not negligible. For stronger 
magnetic fields, i.e. higher Hartmann numbers, 
the behavior is similar. More details can be 
found in [1 7]. 


Design of a MHD mockup experiment 
of a water cooled lead lithium TBM for 
ITER 


Another concept for an ITER TBM is the water 
cooled lead lithium (WCLL) blanket [18], where 
PbLi is employed as breeder and neutron mul- 
tiplier. Pressurized water cools the first wall 
and breeding zone (BZ), and EUROFER steel 
is used as structural material [19]. The design 
of the WCLL TBM consists of 16 breeder units 
(BU) arranged in two columns with a complex 
structure of PbLi and water manifolds (see Fig- 
ure 10). A system of water cooling pipes placed 
inside BUs removes the heat generated in the 
liquid metal. These water pipes are crossing 
the PbLi manifolds, and they occupy a signifi- 
cant fraction of the manifold cross section. In 
addition, stiffening plates inside the PbLi man- 
ifolds lead to repeated reduction of the cross 
section and redirect locally the liquid metal 
flow. The meandering and constricted flow 
paths increase the PbLi velocity and pressure 
drop in the manifolds. 
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Figure 10: Design of the WCLL TBM as provided by CEA 
[18]. View on details such as liquid metal manifolds, water 
pipes, breeder units, baffle plates and stiffening plates. 


Objective of the planned experiments is to 
achieve knowledge about pressure drop in the 
TBM of the WCLL blanket and to determine the 
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distribution of the PbLi flow from the manifold 
into the breeder units. With the aim to investi- 
gate the liquid metal flow in a blanket geometry 
that is most realistic and comparable to the 
original concept foreseen for ITER, a scaled 
MHD mockup has been designed on basis of 
the present WCLL TBM design [18]. In order to 
achieve high Hartmann numbers in a range rel- 
evant for the application, the mockup geometry 
should be as large as possible. Nevertheless, 
the mockup has to fit into the gap of the magnet 
available at the MEKKA facility. To keep the 
geometry as close as possible to the real de- 
sign and as large as possible, one column of 8 
BUs with distributing and collecting manifolds 
is considered. With a scale of 1:2.5, such a 
model geometry fits well into the magnet. It has 
all geometric features and details, which are 
necessary for performing meaningful MHD ex- 
periments. 


For the experiments in the MEKKA laboratory, 
the WCLL mockup will be embedded with hor- 
izontal orientation into the magnet. Figure 11 
shows a transparent view of the mockup for 
visualization of the PbLi flow path showing the 
main dimensions. The liquid metal is fed into 
the mock up and removed from it through cir- 
cular pipes, which are connected to the exist- 
ing loop in the laboratory. The PbLi is distrib- 
uted and collected by two manifolds into and 
out of the BUs. Each BU is fed through a small 
window in the back wall of the manifolds. The 
PbLi flow is redirected at the first wall and 
guided back into the outlet manifold. 
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Figure 11: Transparent view of the WCLL MHD mockup 
design with main dimensions (in mm) and indicated PbLi 
flow paths. 


The WCLL TBM mockup has been simplified 
to meet manufacturing concerns and to reduce 
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fabrication costs, while keeping the geometry 
for the liquid metal flow as foreseen in [18]. The 
test section has been designed with the objec- 
tive to reduce the number of parts and the 
number of welds, respectively. An optimized 
and well-placed number of electron beam 
welds minimizes welding distortions. Costs for 
fabrication remain moderate due to the low 
number of welds. 


Details concerning the water flow inside the 
blanket walls have been omitted since it is ex- 
pected that they will not affect the liquid metal 
flow. Water pipes, representing a partial block- 
age for the liquid metal flow in the manifolds 
are all present in the mockup. 


An exploded view of the mockup in Figure 12 
shows the major parts. The main body consists 
of a single piece, which forms the walls of all 
BUs, including first wall, stiffening plates, back 
plate and walls separating the two inlet and 
outlet manifolds for the PbLi flow. The breeder 
units are closed from both sides and the mani- 
fold is closed by two back plates. The water 
pipes blocking part of the liquid metal domain 
are simulated by dummy parts of solid mate- 
rial. All parts are made of 1.4571 austenitic 
steel that has good compatibility with the used 
model fluid NaK, and does practically not dis- 
turb of the magnetic field. 
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Figure 12: Exploded view of the WCLL MHD mockup with 
major components: main body, top and bottom walls, 
back plates as cover for the manifolds, and dummy parts 
that mimic the water pipes. 


Pressure differences will be measured be- 
tween several points of the mockup using a 
piping system between pressure taps and the 
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pressure transducers. The transducers are lo- 
cated outside of the magnet. The piping sys- 
tem is connected with an array of 30 computer- 
switched valves, which control the measuring 
sequence. To avoid additional pipes for drain- 
ing and venting during the filling and emptying 
process, the positions of the pressure taps 
have been chosen to also satisfy the function 
of feeding, draining and venting. They are lo- 
cated at the lowest and highest positions of 
each breeding unit in the breeding zones and 
manifolds (see Figure 13). 
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Figure 13: Connections with liquid metal loop and location 
of pressure taps. 


The aim of the planned experiment is to assess 
MHD pressure losses when the liquid metal 
passes around obstacles (pipes) and through 
contractions and expansions. Flows in different 
parts of the module are electrically coupled 
since currents induced in one part may leak 
across electrically conducting walls into neigh- 
boring channels. Electromagnetic flow cou- 
pling and flow partitioning among BUs is inves- 
tigated by measurements of the electric 
potential distribution on the walls. 


First experiments in a test section 
related to magneto-convection in 
water-cooled lead lithium blankets 


The numerous cooling tubes immersed in the 
liquid metal of a WCLL blanket module create 


large temperature gradients such that the liq- 
uid metal flow, which is partially obstructed by 
these obstacles, results from the combined in- 
teraction of electromagnetic force, buoyancy, 
and pressure gradient. To improve the under- 
standing of magneto-convective flows, a sim- 
plified mockup is used to study buoyant MHD 
flows in cavities with internal obstacles and to 
provide a first experimental database for MHD 
flows with heat transfer in a WCLL-like geom- 
etry. 


Figure 14: Instrumented test section with transparent 
walls for preliminary water experiments. 


The experimental test section is shown in Fig- 
ure 14. It was designed to be used with GalnSn 
as a model fluid and to fit within the gap of the 
electromagnet present in the MEKKA labora- 
tory. It consists of a rectangular box made of 
PEEK plastic through which two parallel cop- 
per pipes serving as a heat source and a heat 
sink are inserted. Each pipe is connected to its 
own temperature-controlled circuit such that 
they can be kept at constant temperatures T, 
and T, during the experiments. This set-up pro- 
vides the necessary differential heating to es- 
tablish a horizontal temperature gradient that 
drives the buoyant flow. The pipes are coated 
with a 2 um layer of silicon carbide to electri- 
cally insulate the copper from the liquid metal 
and protect them from corrosion. The test sec- 
tion is fully instrumented to record simultane- 
ously local velocities, electric potential at the 
walls and temperature at several locations in 
the fluid. In particular, the temperature distribu- 
tion in the center of the cavity is measured by 
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an 11-point thermocouple probe uniformly dis- 
tributed across the height of the box. Further 
details are given in [20]. 


Before proceeding with the liquid metal exper- 
imental campaign, preliminary experiments 
have been performed with water to test the in- 
strumentation and to obtain results for the hy- 
drodynamic limit where no magnet field is ap- 
plied. Besides, using a transparent medium 
also enables full flow visualization by optical 
techniques. For that purpose, two walls of the 
box were replaced by Plexiglas panels and the 
water was seeded with fluorescent tracers so 
that Particle Image Velocimetry can be applied 
(Figure 15). 


Figure 15: Laser sheet produced by a Nd:YAG laser illu- 
minating the test section during PIV measurements. 


Experiments have been performed at various 
temperature difference AT = T, - T, imposed 
between the pipes. Results obtained show that 
the buoyant flow results in a thermal stratifica- 
tion with the hot fluid staying on the top and the 
cold fluid on the bottom as seen on the dimen- 
sionless temperature profiles measured at the 
central probe and plotted in Figure 16. 
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Figure 16: Non-dimensional temperature profile meas- 
ured at the center of the cavity. 


After finishing the preliminary water campaign, 
the test section will be filled with liquid metal 
and placed in the magnet to perform experi- 
ments for various magnetic fields and temper- 
ature differences. 


Further work 


In addition to the topics described above, the 
MHD group at IKET KIT contributed in the re- 
porting period 2019 with scientific papers to the 
development of fusion technology [21] [22] 
and to fundamental studies of MHD flows [23] 
[24] [25] [26] [27] [28] [29]. 
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3-D Simulation of Fuel Assembly Blockage in MYRRHA 


The paper is based on the presentation by X.N. Chen et al. at TCADS-4, Fourth Interna- 
tional Workshop on Technology and Components of Accelerator-Driven Systems, Ant- 


werp, Belgium, 14-17 October 2019. 


Xue-Nong Chen*, Andrei Rineiski 


Abstract 


In the framework of the KIT and SCK-CEN 
R&D cooperation and as continuation of earlier 
studies performed for the EURATOM FP7 
MAXSIMA project, accidental transients 
caused by a single fuel assembly (FA) block- 
age were simulated with the SIMMER-IV (3-D) 
code for the MYRRHA core, while assuming no 
power variation during the accident. A 7-FA 
model that includes mesh cells for inter-wrap- 
per gaps between FAs was applied, with the 
blockage of the central FA. Sensitivity anal- 
yses on the gap flow rate, fuel chunk jamming 
fraction, insulator pellet material were per- 
formed in order to identify a conservative case 
that maximises the chance of damage propa- 
gation from the blocked FA to the neighbouring 
ones. All calculations with different options and 
parameters did show the same sequence in 
the blocked FA, including melting of pin clad- 
ding, fuel pellet failure, small can-wall break- 
up, steel particle and fuel chunk accumulation 
leading to additional blockages, and large can- 
wall break-up. Finally, fuel chunks are swept 
out from this FA through the inter-wrapper 
gaps. In the calculations performed for several 
tens of seconds and longer, no canwall break- 
up in the neighbouring FAs has been ob- 
served. Nevertheless, different options for sim- 
ulation of the insulator break-up lead to signifi- 
cantly different results in the later phases of 
calculations. If the insulator pellet breaks up 
when the cladding is lost, a fuel/steel blockage 
is formed, which results in a large canwall 
break-up, but this blockage is dissolved as 
soon as the upper steel structure melts. If no 


insulator pellet breaks-up, the fuel/steel block- 
age is kept in place by the ceramic insulator, 
which has a very high melting temperature. 
This observation supports the use of an insu- 
lating material with low melting temperature; 
this option may prevent or reduce blockage of 
fuel/steel particles coming from failed pins that 
eventually may prevent or reduce the risk of 
damage propagation to the intact core. 


Keywords 


Reactor safety study, LBE cooled reactor, 
MYRRHA, SIMMER code, 3-D simulations, 
Fuel assembly blockage accident 


Introduction 


The MYRRHA reactor designed and devel- 
oped by the Belgian Nuclear Research Centre 
(SCK-CEN) was investigated in the European 
7% framework project MAXSIMA [1] and earlier 
projects. The present work is dedicated to 
analyses of fuel assembly (FA) blockage acci- 
dents in MYRRHA studied at KIT in coopera- 
tion with SCK-CEN. 


SIMMER-III/IV is a two-/three-dimensional, 
multi-velocity-field, multi-phase, multicompo- 
nent, Eulerian, fluid-dynamics code system 
coupled with a structure model including fuel 
pins, hexcans, etc., and a space-, time- and 
energy-dependent transport theory neutron dy- 
namics model [2, 3]. The overall fluid-dynamics 
solution algorithm is based on a time-factoriza- 
tion approach, in which intra-cell interfacial 
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area source terms, heat and mass transfers, 
and the momentum exchange functions are 
determined separately from inter-cell fluid con- 
vection. In addition, an analytical equation-of- 
state (EOS) model is available to close and 
complete the fluid-dynamics conservation 
equations. The code was originally developed 
and applied for severe accident analyses of 
fast sodium cooled reactors. However, the phi- 
losophy behind the SIMMER development was 
to generate a versatile and flexible determinis- 
tic tool, applicable for the safety analysis of var- 
ious reactor types with different neutron spec- 
tra and coolants, up to the new accelerator 
driven systems for waste transmutation. In par- 
ticular, SIMMER was extended at KIT to stud- 
ies of heavy-metal-cooled, gas-cooled and 
molten-salt reactors. 


The SIMMER IV code, validated by KIT exper- 
imental KIT [6], see Fig. 1, was applied for the 
described studies. 
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Figure 1: Comparison of the coolant and rod tempera- 
tures from SIMMER-IV simulation and experiment [6]. 


A geometrical model for the chosen 7-FA con- 
figuration of a recent design version of the 
MYRRHA critical reactor has been set up, 
where the interwrapper gaps have been explic- 
itly taken into account. Central FA blockage ac- 
cidents are studied, where the configuration of 
central blocked FA surrounded by other 6 FAs 
are asymmetric. Results are presented and an- 
alysed in this paper. 
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SIMMER-IV modeling for the 7-FA 
conficuration 


Geometrical Modeling 


The coolant flow through gaps (inter-wrapper 
flow) between FAs cools the FA canwall and 
therefore can delay or even prevent the FA 
canwall failure and its propagation. As SIM- 
MER was developed initially for simulations of 
core disruptive accidents (CDAs) in sodium 
fast reactor and in simulations of such acci- 
dents explicit — i.e. via allocation of special 
meshes in plane - treatment of gaps between 
FAs is obviously not important, the the mainly 
used modelling approach with SIMMER was 
not to not consider gaps explicitly. However 
when the SIMMER code is applied to other 
scenarios than CDAs, the modelling of gaps 
may reduce the modelling uncertainties. In the 
past, the gaps were modelled explicitly in a 2D 
case providing quite different simulation results 
as compared to the implicit option [4, 5]. More- 
over the gaps were explicitly modelled also in 
3D with SIMMER-IV for an experimenat case 
[6]. Due use of the XYZ-geomtery in SIMMER- 
IV, the hexcan geometry is transformed into a 
rectangle cubic geometry. Taking 7 FAs as an 
example, the hexcan geometry could be con- 
verted into a rectangular one, as shown in Fig. 
2. Each FA has to be divided into 2 parts with 
suitable parameters. The upper part "(p)" con- 
tains left, right, and back canwall. The lower 
part "(q)" contains left, right, and the front can- 
wall. Additionally the gaps between the subas- 
semblies need to be explicitly modelled as part 
“(r)”. As long as the explicit meshes are set up 
in the computational system, SIMMER code 
treats the flow in the gaps with the hydraulic 
diameter based on the actual geometry, the 
heat transfer between the canwall (steel) and 
the flow (LBE) is calculated with the Nussult 
numbers based on empirical heat transfer cor- 
relations. 


Power Amplitude and Distributions 


The MYRRHA critical core has 108 FAs. 7 FAs 
have been selected for the study. The selected 
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7 FAs are numbered and their power peaking 
factors are indicated in Fig. 3. The selection of 
the central blocked FA is representative for a 
heterogeneous combination of subassemblies, 
with a maximum risk of a local failure of the 
canwall of the blocked FA. In this case the pos- 
sible following fuel sweep-out is anisotropic 
with a high probability of damaging the neigh- 
boring FAs. The FA-wise and axial power dis- 
tributions are imposed on each single assem- 
bly according to the reference results [7]. The 
power is assumed to be constant during the 
blockage transients. The thermal boundary 
conditions outside the system are adiabatic. 


KA 


Figure 2: Explicit meshes for gaps flow between subas- 
semblies using SIMMER-IV code. 


Figure 3: 7 fuel assembly arrangement with power peak- 
ing factors 


The steady states at the BOC and EOC condi- 
tions are quite similar, where only the peaking 
factors at BOC are slightly higher than those at 
EOC, mainly due to different axial positions of 
CRs. Therefore, to be conservative, the BOC 
condition was used for steady state and transi- 
ent analyses. 


Key steady-state thermohydraulic parameters 
calculated by SIMMER-IV are calculated for 
every FA, which match the design values [8]. 
The coolant inlet temperature is 270 °C. The 


average thermal power per FA is 0.9259 MW. 
The pressure boundary conditions are set up 
and the orifice coefficients for the seven FAs 
and gaps (in the FA upper position) are ad- 
justed, so that certain coolant flow rates (veloc- 
ities) are achieved in FAs and gaps respec- 
tively. It is noticed that the flows in the gaps are 
heated up by the can wall. 


FA Blockage Transients 


FA blockage conditions and parametric 
studies 


The FA blockage is a postulated envelope sce- 
nario that evolves from local pin failures which 
propagate at a FA scale. The blockage is sim- 
ulated by a reduction of the FA inlet flow rate. 
Its flow rate is reduced from 72 kg/s to 3.2 kg/s 
(i.e. 95.6% blockage) for FA1, the central FA, 
as itis more realistic to assume a porous block- 
age, allowing a residual cooling flow. The fuel 
particle diameter is 5.6 mm. For these calcula- 
tions the fission gas blowout model is deac- 
tivated, because we assume, the fission gases 
contained in the pin gap are all released before 
the blockage extends to the whole FA flow 
area. 


Different parameters, such as gap velocity, in- 
sulator material, particle packing fraction limit 
in the particle jamming model, and other break- 
up model options, were employed, some re- 
sults are given in the following. The reduction 
of the gap velocity has no signicant effects on 
the clad melting, but it causes a numerical in- 
stability problem. The simulations with different 
insulator materials as steel and B4C give very 
similar results. The increase of the particle 
packing fraction limit shows no significant dif- 
ferences in the calculated results except a 
slightly earlier canwall breakup and larger par- 
ticle blockage coefficients. 


19 


Group: Transmutation 


Time=32.400963 [38609 731.96 1871.81 
E Z 
< < 
543.006 543.012 
3 
Radial 
TLK3 TLK3 
Time-20.000734 Time-34 000404 [41948 
7 Gas Gas 7 Gas 
Steel Steel Steel 
Pellet Pellet Pellet 
Li Ll Ll 
L2 L2 L2 
L3 L3 3 
L4 I4 L4 
L5 LS Ls 
= L6 L6 3 L6 
E L7 L7 z L7 
Žž < ^ 
B Crust Crust J Crust 
Cavity Cavity Cavity 
Other Other Other 
3 3 3 
Radial Radial Radial 
2104.5 Time-42. 800018 [1188158 IJ (K=4) 2319.45 Time-47 816017 [2103206 2278.67 
3 E z 
E Z E 
2 2 2 
543.02 543.021 543.019 
Gas Gas Gas 
Steel Steel Steel 
Pellet Pellet Pellet 
LI Ll L1 
L2 L2 L2 
3 L3 L3 
L4 L4 L4 
L5 L5 LS 
L6 = L6 3 L6 
L7 E] L7 E L7 
Crust Crust S Crust 
Cavity Cavity Cavity 
Other Other Other 


Radial Radial Radial 


Figure 4: Fuel pin breaks up and fuel chunk blockage: coolant temperature TLK3 and material distribution, where 
L1 liquid fuel, L2 liquid steel, L3 liquid sodium, L4 fuel particle, L5, steel particle, L6 control particle, L7 fuel chunks. 


20 


Group: Transmutation 


FA blockage transient results in a typical 
case 


In this paper we just present the calculated 
blockage results in a typical case, i.e. gap ve- 
locity = 0.22 m/s, insulator material is B4C, the 
particle packing fraction limit APJ = 0.8. 


Fig. 4 shows the fuel pin breaks up, coolant 
temperature (TLK3) and material distributions 
after 20 seconds, the starting time of the block- 
age. As long as the blockage occurs, the cool- 
ant temperature reaches about 1000 °C after 
less than 10 s. As the cladding temperature 
reaches the steel melting point limit (1430 °C), 
the cladding starts to melt. The fuel dispersion 
starts when the cladding is broken up after 12.4 
seconds of instantaneous blockage. A large 
fuel chunk blockage takes place then aflter 14 
seconds at the upper part of the active core. 
Then the canwall starts to break up after 16.2 
seconds. The fuel chunk blockage is not dis- 
solved during the transient calculated. The 
canwall breakup takes place only at the middle 
small FA cells. Afterwards, the fuel in this FA 
escapes from this break-up through the gap to 
the upper part of the core. No neighboring can- 
wall break-up takes place, i.e. the neighboring 
FAs stay intact during the transient. Neverthe- 
less, additional blockages in the interwrapper 
space cannot be excluded, as they cannot be 
detected with the current model. It is noticed 
that, if the blockage is a little bit smaller, the 
canwall in the blocked FA would not break up, 
thanks to the cooling provided by the flow in 
gaps, as reported in [9]. This means that the 
gap flow can prevent the canwall break-up and 
its break-up propagation. 


Conclusions 


Calculations with different options and param- 
eters show the same sequences, i.e. the pin 
melting, the small canwall break-up, fuel/steel 
particle accumulation and its blockage and the 
canwall break-up in the blocked FA and finally 
fuel particle swept out from the blockage FA 
through the inter-wrapper gaps. The conserva- 
tive calculations show there is no canwall 


break-up in the neighboring FAs, even if addi- 
tional studies need to be carried out to exclude 
the formation of fuel blockages in the inter- 
wrapper space. The variation of the maximum 
packing fraction APJ shows that it has only 
slight effects on the later sequences of particle 
blockage, which does not change major results 
and conclusion. The variation of the interwrap- 
per gap velocity shows that it has some effects 
on the clad melting time and later sequences. 
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JIMEC experiments to investigate jet impingement on 
a core catcher bottom and LIVE2D 2-Layer experiment 


Xiaoyang Gaus-Liu, Thomas Cron, Beatrix Fluhrer, Rene Stángle, Mike Vervoortz, Thomas Wenz 


Introduction 


In 2019, the research activities in IKET-SAR 
team were focused on several European pro- 
jects. Some final work had to be done for the 
European SAFEST project, coordinated by 
IKET-SAR team, which ended officially at the 
end of 2018. 


The main activities in 2019 were concentrated 
on the European ESFR-SMART project. Within 
this project two large-scale JIMEC experi- 
ments have been performed to investigate the 
thermal ablation kinetics of an internal core 
catcher material in a SFR reactor. Besides this, 
the planning work for LIVE-ESFR tests to study 
the interaction between the corium simulant 
and the sacrificial simulant of the core catcher 
started. It has been decided to construct and 
build a new test vessel with down-scaled ge- 
ometries similar to SFR core catcher design. 


Another activity is the finalization of LIVE2D 2- 
Layer experimental and analytical analysis, 
which is a main experimental task in the H2020 
IVMR project, which ended in November 2019. 


JIMEC experiments to investigate 
jet impingement on a core catcher 
bottom and ablation process 


Background and Objectives 


The actual safety design of a Sodium Fast Re- 
actor (SFR) in the case of a postulated severe 
accident incorporates to remove the corium 
from the core by corium transfer tubes and to 
collect it in the lower head in an in-vessel core 
catcher. It is assumed that at first a metallic co- 
rium melt jet would impinge on the core catcher 


surface and could ablate the core catcher ma- 
terial. Experimental data is needed to simulate 
this ablation behaviour of a long duration melt 
jet impinging the core catcher material. A par- 
ticular behaviour can be studied when a molten 
pool is created (“pool effect”) at the impact 
point that could reduce the heat transfer at the 
jet — material interface. This phenomenon has 
been studied very little in the past. Therefore, 
the IKET-SAR team has adapted the existing 
MOCKA test facility to perform two JIMEC (Jet 
Impingement on Metallic Core Catcher) exper- 
iments in the frame of the European ESFR- 
SMART project. JIMEC-1 and JIMEC-2 tests 
investigate the characteristics of the interaction 
of a metallic melt jet with the core catcher bot- 
tom plate material in a SFR reactor design with 
prototypical material. The objectives of the ex- 
periments are to deliver experimental data on 
the interaction of melt jet parameter and ero- 
sion dynamics. The melt jet parameters were 
jet temperature, jet velocity and jet diameter. 
The erosion dynamics in the core catcher bot- 
tom is the erosion velocity and the timing of 
pool effect. The experimental results will be 
used for developing new correlations which 
could be used in codes for simulation of the ab- 
lation kinetics for SFR core catcher concepts. 


Experiments 


The two JIMEC experiments have been per- 
formed in summer 2019 in the adapted 
MOCKA test facility, Figure 1. To simulate the 
core catcher bottom a test substrate was used 
composed of a oylindrical formed stainless 
steel block (1.4301) with 416 mm thickness 
and 425 mm in diameter. A matrix of thermo- 
couples (TCs) was implemented into the test 
substrate to record the erosion by the melt jet. 
About 1000 kg of metallic melt was produced 
by thermite reaction in the separate reaction 
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crucible above the test substrate. The metallic 
melt is designed to have the same composition 
as the metallic test substrate. The outflow 
opening at the bottom of the crucible was 40 
mm in diameter in JIMEC-1 and 30 mm in 
JIMEC-2. The thermite reaction produced also 
a lighter oxide melt consisting mainly of Al2Os, 
Cr2Os and CaO. This melt was redirected by a 
pouring spout to an oxide melt collector. 
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Figure 1: Picture of JIMEC-2 test set-up in top picture 
and the thermocouple instrumentation in the substrate 
in the bottom sketch. 
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For the jet outlet diameter of 40 mm in JIMEC- 
1, the duration of the metallic melt jet was 
about 31 s from start of outflow to the end of 
metallic melt jet. In JIMEC-2 the duration was 
55 s due to the smaller jet outlet diameter of 30 
mm. In JIMEC-1, the start of pool effect was 
about 18 s after the first melt reached the sub- 
strate indicated by the stop of splashing of the 
melt, Figure 2. In JIMEC-2, the start of pool ef- 
fect was detected about 15 s after the first melt 
reached the substrate surface. Infrared pic- 
tures of JIMEC-2 before and after start of pool 
effect are shown in Figure 3. To protect the en- 
vironmental apparatus in against the hot 
splashing melt in JIMEC-2 test, a half-cylinder 
concrete pipe was positioned upon the sub- 
strate with an opening at the jet flow position. 


Figure 2: Video pictures of melt jet before and after start of 
pool effect in JIMEC-1 


Figure 3: Infrared pictures of JIMEC-2 before and after 
start of pool effect. 
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The shape of the jet remained in both tests 
most time coherent, which means there was no 
breakup of the jet flow. The jet outflow diame- 
ter enlarged during the jet release process 
since the outlet nozzle made of ZrO2 was grad- 
ually ablated by the metallic and also by the 
oxide melt. The post-test diameter is 54 — 58 
mm for JIMEC-1 and about 46 mm for JIMEC- 
2. Analysis from high-speed video of JIMEC-1 
shows that before the pool effect the melt jet 
velocity is 4.3 - 4.6 m/s and after the pool effect 
about 5.0 m/s. The theoretical velocity is about 
4.4 m/s at the outlet nozzle. The large velocity 
at the later phase could be an effect of the en- 
larged outlet diameter. The metallic melt jet 
temperature in both tests was in the range of 
2000 - 2100 °C measured by a pyrometer. The 
ablation velocity in the substrate can be ob- 
tained based on the thermocouple signals. Fig- 
ure 4 shows the ablation depth in the substrate 
versus time for JIMEC-1 and JIMEC-2. 


Comparison of ablation depth in substrate 
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Figure 4: Comparison of ablation depths for JIMEC-1 and 
JIMEC-2 


The ablation rates within 30 mm diameter in 
the substrate are almost identical for both ex- 
periments. The ablation velocity in the center 
of the substrate is about 17 - 18 mm/s before 
start of pool effect and decreases to about 7 
mm/s after start of pool effect. An influence of 
the different jet diameter on the ablation be- 
haviour is then detected in the outer regions of 
the test substrate. For JIMEC-1 with 40 mm jet 
diameter at the beginning, the radial ablation 
proceeds faster than for JIMEC-2 with 30 mm 
jet diameter. 


According to the ablation velocity, the contour 
of ablated pit vs. time can be roughly esti- 
mated, as shown in Figure 5. Due to the lacks 
on measuring position at the outer bottom part, 
the contour up from 23 seconds can only be 
partially illustrated. The forms of the contours 
indicates the initiation of pool effect in17 sec, 
leading to a slowdown of the axial ablation, 
however a high radial erosion rate at the lower 
part, and thus this led to a gradual transfor- 
mation of a conic pit to a cylindrical pit. 
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Figure 5: Progression of ablation contour in the JIMEC-1 
substrate 


The test substrates for JIMEC-1 and JIMEC-2 
have been cut after the tests, as shown in Fig- 
ure 6. The test substrates of JIMEC-1 and 
JIMEC-2 are molten through by the metallic 
melt jet. In JIMEC-1 a hole was formed in the 
lateral wall above a height of -132 mm from 
the bottom. Up to the same height, a solidified 
metal melt has built in the cavity. The mass of 
both substrates before the tests was about 460 
kg. The remaining mass of the test substrate of 
JIMEC-1 is ~ 240 kg including the solidified 
metal melt pool. For JIMEC-2, the remaining 
mass is about 214 kg. Therefore, about 246 kg 
of the substrate has been molten in IMEC-2 
substrate. 
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Figure 6: Cut of the test substrates of JIMEC-1 and 
JIMEC-2 


Analyses of the LIVE-2D two-layer test 
series on the heat flux focusing effect 
and the thermal-hydraulic character in 
the upper melt layer 


Two series of LIVE2D two-layer tests with 3 up- 
per layer thicknesses and different surface 
boundary conditions were carried out in 2017. 
The experiment reveals major thermos-hy- 
draulic characteristics of the upper light melt 
layer during the transient and steady states as 
well as the strong dependence of the heat flux 
focusing effect on the upper boundary cooling 
condition. The final analysis of the two test se- 
ries comprises the study the heat flux focusing 
effect upon different boundary conditions and 
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the characterization of the thermal-hydraulic 
feature of the upper layer. 


Strong heat focusing effect was observed at 
the wall in the upper melt layer in LIVE2D-SO1 
test, which had a hot upper atmosphere. And 
no heat flux focusing effect appeared when the 
melt upper surface was rigidly cooled, which 
was the upper boundary condition of the sec- 
ond test series (LIVE2D-SOTC). Figure 7. 
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Figure 7: Heat flux at the vessel wall with 110 mm upper 
layer. Top: during SO1 test; bottom: during SOTC test. 


The experiment results demonstrate further 
that without effective cooling at the upper 
boundary, the melt circulates globally from the 
hot central region radially toward the cooled 
wall boundary, resulting a large diversion of 
bottom boundary temperature, as shown in 
Figure 8. Whereas dimensionless temperature 
0' = (T —Tmin)/(Tmax — Tmin), Z: distance to 
the bottom and L: the layer thickness. Upon 
this turbulent flow character, the generally ap- 
plied Globe&Dropkin correlation, describing 
convective heat transfer upon a uniform hot 
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bottom plate to the cold bulk melt in vertical di- 
rection, is not suitable for the upper layer heat 
transfer with limited upper heat transfer rate. 
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Figure 8: Dimensionless temperature vs. vertical position of a 110 mm upper layer inSO1 test 


In addition, the specific boundary combination 
of the upper layer without strong top cooling 
doesn't exactly correspond the sidewall heat 
transfer situation of the widely applied Church- 
ill&Chu correlation, which describes the heat 
transfer of a heated vertical wall to the cold en- 
vironment. In Figure 9 the LIVE2D experi- 
mental results are plotted in comparison with 
the two Churchill&Chu correlations, and the 
experimental results show a considerably 
lower convective heat transfer capability than 
the Churchill&Chu correlations. 
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Figure 9: Comparison of Nu at the heat transfer at the upper layer sidewall 
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Detailed investigations on flow boiling of water up to the 


criticalheat flux 


Stephan Gabriel, Giancarlo Albrecht 


Introduction 


Multiphase flows occur in various fields of pro- 
cess engineering and energy technologies. 
Highly complex phenomena such as flow boil- 
ing up to the critical heat flux (CHF) or counter 
current flow limitation CCFL arise from a whole 
series of phenomena such as convection, bub- 
ble formation, evaporation, condensation, 
drop/bubble fragmentation or coalescence and 
many more. In order to to understand the be- 
haviour of two-phase flows in real technical 
systems in detail, it is necessary to consider 
both the individual phenomena and the more 
complex phenomena in the experiment and to 
advance the development of simulation pro- 
grams with detailed data sets. 


The Group Multiphase Flows (MPF) works 
therefore on investigations at different scales 
from experiments on single effects [4] to exper- 
iments on complex phenomena like CHF. 
Since the required measurement technology is 
in many cases not commercially available, the 
development of suitable sensors, for example 
for the acquisition of detailed phase distribu- 
tions, is also part of the task. 


Results 


In 2019, the focus of work was on the construc- 
tion of the high-pressure test facility COSMOS- 
H on the one hand, and on the other the com- 
pletion of measurements on a rod bundle test 
section of the low-pressure test facility COS- 
MOS-L. As a further important work package, 
extensive work was carried out on modernizing 
and upgrading the laboratory infrastructure in 
view of the upcoming commissioning of COS- 
MOS-H. This includes the procurement or re- 


commissioning of laboratory equipment such 
as tensiometers, microscopy and measuring 
instruments for determining water quality for 
the tests and as a precondition for the tests the 
upgrading of the crane technology for the new 
test track. In addition, the dismantling of the 
meanwhile 20-year-old test facility DISCO was 
carried out to create urgently needed space for 
COSMOS-H. 


COSMOS-L 


At the COSMOS-L test facility in 2019, further 
measurements were made to complete a da- 
taset on a rod bundle and then a new test track 
with a flat vertical heater was developed. A 
measurement data set for the rod bundle con- 
sisting of 5 Zircalloy tubes, which was created 
within the framework of the joint project 
NUBEKS [2,3], was supplemented by further 
measurement points and measured variables. 
The measurements characterize the behaviour 
starting from flow boiling up to the boiling crisis 
(CHF). As can be seen in Figure 1, the test 
section has a rod bundle consisting of five sep- 
arately heatable cladding tubes. During the 
tests, either all tubes were heated in order to 
obtain the most realistic phase distribution 
around the central tube, or only the central tube 
was heated in order to concentrate the instru- 
mentation there. The results presented below 
were obtained by experiments with only one 
heated tube. 


The test results presented in the following refer 
to a static pressure of 1200 mbar, a mass flow 
density of 50 kg/(m?s) and a test track inlet 
temperature of 80?C. This corresponds to an 
inlet subcooling of approx. 25°C. As shown in 
Figure 2, bubble boiling already occurs at a 
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heating power of 4 kW, which becomes more 
and more intensive with increasing heat flux. 
The critical heat flux is reached at this bound- 
ary conditions at a heating power of approxi- 
mately 10.8 kW. 
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Figure 1: Experimental setup of the bundle measure- 
ments 


Analogous to the previous experiments 2018, 
the experiments were repeated several times, 
here a total of 29 times, and a frequency distri- 
bution of the CHF value was generated (see 
Figure 3). The measured values related to the 
heated surface thus result in a critical power of 
CHF mean 21.1257 MW/m? with a standard devi- 
ation of o 20.013 MW/m?. 
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Figure 2: Flow pattern as a function of heat output 


In addition, the locations where the CHF oc- 
curred on the cladding tube could be calcu- 
lated by correlation of the three thermocouples 
closest to the hotspot [5]. For this purpose, the 
heated pipe in the test section is divided into 
triangles between each three thermocouples. 
The measured temperature values show differ- 
ent signals depending on how close they are to 
the location of the CHF. Characteristic here are 
on the one hand the waiting time, i.e. the run- 
ning time of the signal, and on the other hand 
also the temperature gradient. 


After the basic function of the principle was 
demonstrated on a sheet metal of 40 x 40 cm? 
with a propane gas burner as heat source, the 
CHF events of the COSMOS-L measurements 
were evaluated concerning the frequency dis- 
tribution shown above. It is evident that the 
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Figure 3: Frequency distribution of the critical heat flux after 29 tests 
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method has room for improvement in terms of 
precision. However, it can certainly be used to 
determine during a test whether the CHF 
events always occur at the same location, 
which would be an indication of damage to the 
heating tube. 


As expected, most CHF events could be iden- 
tified at the upper end of the 326 mm long 
heated tube. It was found that although there 
are accumulations of CHF events but they 
rarely occur directly one after the other [5]. 


Following the rod bundle experiments, the de- 
sign and construction of a new test arrange- 
ment for experiments within the project KEK- 
SIMA was started, in which, together with the 
Institute of Applied Thermofluidics (IATF), the 
effects of transient mass flows on the critical 
heat flow density will be investigated. For this 
purpose, a test section with a plate-shaped 
heating element was constructed in 2019 (see 
Figure 5). In addition to a powerful heater, the 
test section also has three window elements 
for observing the flow and numerous sensors 
for recording all relevant boundary conditions. 
A precondition for the upcoming tests is an in- 
crease in performance and an extension of 
COSMOS-L by numerous pressure, tempera- 
ture and mass flow sensors. In addition, a by- 
pass to the test track will be installed in the 
coming year in order to ensure stable system 
operation and thus reproducible test boundary 
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conditions even with fluctuating mass flows in 
the test track. 
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Figure 5: New test section concept with flat heater to in- 
vestigate the critical heat flow density with flucuating 
mass flows 


COSMOS-H 


The thermohydraulic test facility COSMOS-H 
is a high pressure high temperature water loop 
being built for the investigation of boiling phe- 
nomena and other flow phenomena that can 
occur in thermal power plants. In contrast to 
COSMOS-L, the loop will achieve non-scaled 
test conditions comparable to a boiling water 
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Figure 4: Local distribution of CHF events on the middle pipe [5]. 
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reactor or pressurized water reactor (up to 170 
bar, 360°C) during test operation. The work on 
the construction of the COSMOS-H test facility 
was continued as planned in 2019. Numerous 
subsystems were completed and some have 
already been commissioned (see Figure 6). 
These include, for example: 


e The power supply for the control system 
and all electrical consumers was com- 
pleted 


e The two cooling loops of the facility were 
completed (cooling capacity 1.8 MW) 


e The compressed air supply to drive the 
automatic valves is completed, including 
the emergency supply system. 


e The hardware of the digital control system 
(SPS) as well as the technology of the 
control station is completed. 


e The pumps and valves of the plant are 
now functional. 


e Thecrane technology for the assembly 
and installation of the test track, consist- 
ing of assembly crane, transport trolley, 
crane traverse and hall crane, has been 
completed. 


Further work is now concentrating on complet- 
ing the high-pressure loop and the necessary 
safety systems. 
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Figure 6: Images of various subsystems of the plant un- 
der construction 
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Weizel 


Introduction 


Liquid metals are investigated as heat transfer 
fluids in concentrating solar power systems 
due their excellent heat transfer properties 
[1,2]. In the SOMMER test facility at KALLA, 
lead bismuth eutectic (LBE) is used as heat 
transfer fluid in a so-called thermal receiver in 
the focal point of the solar furnace [3]. During 
2019, the SOMMER test facility has been put 
up into operation and tests were performed un- 
der direct concentrated sunlight conditions. 
The thermal receiver — efficiently cooled by 
LBE — withstood maximum heat flux densities 
of 4 MW/m? without any failure. 


Thermal receiver tests in the 
SOMMER facility 


At the location of the SOMMER facility, peak 
DNI values of up to 950 W/m? have been 
measured. The sunlight is reflected by the he- 
liostat mirror (32 m?) and onto the parabolic 
mirror (16 m?) inside the SOMMER laboratory, 
where the solar power is concentrated. The so- 
lar furnace arrangement without the liquid 
metal loop is shown in Fig. 1. The installed 
shutter blinds in front of the parabolic mirror 
can be adjusted for tests under part-load con- 
ditions. For safety reasons, a curtain has been 
installed that drops into the light passage 
driven by gravity in case of a power outage or 
other irregular conditions. Additionally, the he- 
liostat mirror is moved automatically into stow 
position if there is a power outage. 


Figure 1: Solar furnace arrangement consisting of a helio- 
stat mirror (front, outside) and a parabolic mirror (back, 
inside the laboratory) (Foto: J. Flesch). 


The liquid metal LBE can be pre-heated in an 
electric heater before entering the thermal re- 
ceiver to allow for a wide measuring matrix. Af- 
ter leaving the thermal receiver, the heated liq- 
uid metal is again cooled down before entering 
the pump tank. The pump is a submerged me- 
chanical gear pump delivering 0.1 L/s at 
1000 turns/min and is designed to operate up 
to a temperature of 380 ?C. It is located at the 
highest point of the SOMMER liquid metal 
loop. During down time or in case of an emer- 
gency, the fluid flows down into the sump tank 
due to gravity. The sump tank has capacity for 
130 L of lead-bismuth (Fig. 2). 
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Figure 2: The sump tank of the SOMMER test facility; 
top: insulated (Foto: Amadeus Bramsiepe/KIT); bottom: 
before insulation with electric trace heating (Foto: F. 
Fellmoser). 


The thermal receiver, which is cooled by the 
LBE flow, is installed in the focal point of the 
parabolic mirror. It is designed as a spiral tube 
of 10 mm outer diameter with a wall thickness 
of 0.5 mm. The area being heated by the sun 
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is 100 mm x 100 mm. The inlet and outlet tem- 
peratures of the receiver are measured with 
thermocouples to determine the thermal power 
that is absorbed by the liquid metal flow. In to- 
tal, the thermal receiver has a thermal power 
input of max. 13 kW. Figure 3 shows the ther- 
mal receiver connected to the liquid metal loop. 
The tubes are coated with Pyromark 2500 to 
improve the absorption of the sunlight. A large 
copper plate is installed around the receiver to 
prevent damages on other components in case 
of a tracking error of the heliostat. The copper 
plate will then absorb the heat and the light 
source can be cut off in this time. 


Figure 3: The thermal receiver of the SOMMER test facil- 
ity with copper shield (Foto: Amadeus Bramsiepe/KIT) 


In order to determine the incoming solar power 
on the receiver a heat flux measurement de- 
vice was developed [4—6]. It uses a heat flux 
micro sensor that moves across the focal area 
on a circular path superimposed with a linear 
motion. In the SOMMER test facility, the liquid 
metal loop's position — and thus the thermal re- 
ceiver's position — is fixed. The parabolic mir- 
ror, however, can move along a distance of 
0.5 m driven by a linear motor and thus, the fo- 
cal area can be shifted for the heat flux meas- 
urement. Heat flux density (q’) values of up to 
4 MW/m? were measured and successfully 
cooled by the lead-bismuth flow, as shown in 
Fig. 4. 
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Heat Flux Density Distribution 10° 


Figure 4: Heat flux distribution measured 28.06.2019 with 
open shutter blinds 


In addition to the experimental investigation in 
the SOMMER facility, thermal energy storage 
options for liquid metal were compared [7] and 
a dual-media thermocline storage system was 
theoretically investigated in detail [8]. High- 
temperature thermal energy storage with liquid 
metals was also investigated regarding the ap- 
plication in energy-intensive industries with 
waste heat at high temperatures [9, 10]. Exper- 
iments for such high temperature thermal stor- 
age and liquid metal based process technology 
are ongoing in KALLA [11]. 
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Introduction 


Support of decision making is a form of data 
transformation by analysing complex circum- 
stances and processing large amounts of data 
with the goal to present helpful information to 
decision makers in a simplified and better un- 
derstandable way. Within the terms of this 
broad definition, many methods, algorithms, 
and visualisation means qualify as decision 
support. However, there is no best method in a 
general sense, as all these methods have their 
pros and cons in respect to the specific situa- 
tion they are applied to. In addition, different 
methods can be combined to improve the over- 
all performance of decision support. 


Originally, decision support methods used de- 
terministic data and thus their output was de- 
terministic as well. Consequently, such infor- 
mation leads to binary thinking and hard 
decision making, e.g. if a countermeasure 
strategy is to be chosen, deterministic decision 
support forces the view that one specific strat- 
egy is unconditionally the best and superior to 
all other available strategies. However, de- 
pending on the method used for decision sup- 
port, minor changes in the input data can have 
a huge impact and can cause a different out- 
come. As input data is in general not determin- 
istic but affected by uncertainties, this may re- 
sult in recommendations of suboptimal 
strategies in emergency management. There- 
fore, it is necessary to consider uncertainties in 
decision making methods to improve the over- 
all decision support. 


Uncertainties and Combining 


Uncertainty influencing decision 
making 


Many forms of uncertainty can be identified 
that have an influence on the assessment of an 
emergency and its development over time. The 
following list is certainly incomplete but gives 
an impression of uncertainty types: stochastic 
uncertainties in form of physical randomness, 
epistemological uncertainties by lack of scien- 
tific knowledge, endpoint uncertainties when 
the desired goal endpoint is ill-defined, judge- 
mental uncertainties by defining personal pref- 
erences as facts, and computational uncertain- 
ties by e.g. inaccuracy through numerical 
instability or modelling errors as models are al- 
ways a simplification of the real world and 
therefore limited in one way or the other [1]. 
Sometimes the errors introduced by these un- 
certainties may be small, but as they add up 
they could lead to choosing inferior strategies 
in the end. 


The CONFIDENCE project 


The European project CONFIDENCE (2017- 
2019) aimed to analyse uncertainties and to 
improve the support for emergency manage- 
ment, especially focussing on nuclear acci- 
dents [2]. CONFIDENCE investigated the influ- 
ence of uncertainties on the different phases of 
the full chain of managing a nuclear accident 
beginning with the assessment of data 
(weather, source term), continuing with simu- 
lation of the situation development (dispersion, 
food chain), over analysis of possible counter- 
measure strategies (decision support) up to 
the communication of situation development 
and strategies to the public (social science). 
Within this project the Accident Management 
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Systems (AMS) group of ITES had a leading 
role as coordinator of the project and as leader 
of work package 6, which investigated the in- 
fluence of and coping with uncertainties in de- 
cision making. The work package especially 
focussed on enhancing the existing Multi Crite- 
ria Decision Support (MCDA) tool to handle un- 
certainties as well as Agent Based Modelling 
(ABM) to analyse and better understand the ef- 
fects of uncertainties on the decision making 
process. A special issue of the Radioprotection 
Journal is dedicated to the CONFIDENCE pro- 
ject and will be available in the second half of 
2020. 


Multi Criteria Decision Analysis as 
decision support 


Multi Criteria Decision Analysis (MCDA) co- 
vers various decision support methods, that in 
general provide a ranking on a set of alterna- 
tives by integrating (contradictory) decision cri- 
teria of different scale according to given (per- 
sonal) preferences [3]. The ranking helps de- 
cision makers to choose the best suited alter- 
native, which is frequently, but not necessarily, 
the highest ranked alternative. Considering 
emergency management, MCDA syste-mati- 
cally combines the pros and cons of feasible 
actions to be aggregated into a single numeric 
value, which makes them easily comparable 
between each other. The higher the value, the 
better the according action is rated. In manag- 
ing nuclear emergencies, the actions are in 
general a set of countermeasure strategies like 
"Evacuate people and clean surfaces before 
they return". 


Each action is assigned a ranking value A:,:-:, 
An. The criteria are either quantitative values 
like "Estimated dose" or qualitative values like 
"Public acceptance". The criteria values Cı,---, 
Cm are either simply measured or determined. 
For numerical evaluation, qualitative value 
ranges like {“low’,”high”} have to be mapped to 
quantitative value ranges like {1,2}. Since the 
criteria are typically of different units and 
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scales, the criteria values have to be normal- 
ised onto a unified scale before combining 
them. For this purpose, normalisation functions 
N1,---, Nm, such as e.g. min-max normalisation, 
have to be defined for every single criterion. 
The personal preferences of each criterion are 
represented through weights. The relative im- 
portance of a criterion is reflected in a specific 
normalised weight w1,---,Wm. The normalised 
values of criteria are aggregated in a ranking 
value by using an aggregation method accord- 
ing to their weight. One of the most popular ag- 
gregation methods is the computation of the 
weighted sum, which for each alternative re- 
quires the following computation: 


m 
A, =). We (Gui) Vien 
k=1 


The actions are sorted according to their rank- 
ing values, indicating their order of recommen- 
dation. Ranking and results can be presented 
in multiple ways like e.g. charts, graphs, textual 
report, and others depending on the specific 
requirements of the decision makers (e.g. Fig- 
ure 2, Figure 3). 


Though from the mathematical point of view 
the method is plain and straightforward, there 
are several catches to consider when applying 
it. Firstly, the actions to rank are not generated 
but determined externally, either by another 
tool or by the decision makers themselves. 
Secondly, the determination of feasible criteria 
is also up to the decision makers. As emer- 
gency management is in general subject to a 
group of decision makers respectively adviso- 
ries they have to agree on such a set of criteria 
as a group, finding a common consensus. The 
same holds for determination of qualitative cri- 
teria values, which may be based on personal 
assessment. Finally, the weights are depend- 
ent on personal preference and therefore need 
to be agreed upon within the group. This leads 
to a time consuming setup and intense discus- 
sions among the group members, making the 
MCDA method preferable in situations where 
time is available, e.g. in preparation, training or 
long term recovery decision making. On the 


Group: Accident Management Systems 


other hand, as a benefit, the intense discus- 
sions result in transparency and documenta- 
tion how the ranking and therefore the decision 
was justified. 


Considering uncertainties in MCDA 


The MCDA as described above will process 
deterministic parameters, yet most, if not all, 
scenarios of decision making are affected by 
uncertainties, which requires processing of 
probabilistic parameters. The following section 
describes how this limitation can be overcome. 


Two obvious parts of MCDA can be affected by 
uncertainty: the criteria values and the criteria 
weight values. Such uncertain values can be 
described probabilistically: either as functions 
or histograms. Histograms can be easily achie- 
ved by binning and counting according values, 
e.g. for a histogram of a criterion weight let all 
decision makers provide an integer weight 
value between 1 and 10 according to their pref- 
erence and accumulate the values. On the 
other hand, determining distribution functions 
for criteria values like the “Estimated dose” is 
rather difficult. Yet the important part is not to 
achieve higher accuracy, but to introduce the 
potential variety of values into the ranking, thus 
sensitising decision makers to rather look for 
the most robust solution in all circumstances 
instead of the best solution for one specific 
case. 


MCDA cannot process distribution functions or 
histograms as input values. For this reason, 
ensemble evaluation is applied to overcome 
this limitation. Simply put, from the probabilistic 
MCDA a number of deterministic MCDA are 
generated and evaluated one by one. The de- 
terministic results are combined back into one 
probabilistic result. Because of its simplicity se- 
veral thousand MCDA can be generated and 
evaluated within a second, allowing for large 
sample sets. 


The MCDA tool has been enhanced in that way 
to define probabilistic input, to perform ensem- 


ble evaluation, and to present probabilistic re- 
sults [4]. Figure 1, Figure 2, and Figure 3 give 
an impression on the implemented enhance- 
ments. 
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Figure 1: Definition of a criterion value as normally distrib- 
uted. 
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Figure 2: Overview on probabilistic ranking of 3 alterna- 
tives. In this example "High waste" was ranked first place 
in 53.9 percent of all evaluations of the ensemble set. 
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ternative was ranked better than another one. 
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Agent Based Modelling of the decision 
making process 


The decision making process of stakeholders 
is much more complicated than one would ex- 
pect at first glance, as not only plain numbers 
are considered but also personal preferences, 
experience, prognosis of situation develop- 
ment and behaviour of affected people, etc. 
The decision makers involved in emergency 
management may have different backgrounds 
and may belong to different organisations, 
which is reflected in their individual personality 
and their assessment of the emergency situa- 
tion, thus introducing uncertainties in the deci- 
sion making process. 


To consolidate their different personal deci- 
sions into a single one they have to work to- 
gether, share thoughts, negotiate, and finally 
find a compromise that is acceptable for every- 
one. The best way to understand this process 
would be to interview and observe the decision 
makers while they are confronted with a large 
number of different scenarios. However, this 
would take a considerable amount of time, es- 
pecially of the stakeholder's time, and there- 
fore is not feasible in praxis. A practicable ap- 
proach to address this is to model the be- 
haviour of decision makers and the process of 
decision making. Such a model allows for sim- 
ulation and analysis of a large number of differ- 
ent scenarios, leading to a better understand- 
ing of the underlying uncertainties. 


Agent Based Modelling (ABM) is a program- 
ming paradigm that allows for simulation of (in- 
telligent) individuals, so called agents, and the 
complex interactions between them. More spe- 
cific, a software agent is defined as a computer 
system (program) that is situated in some en- 
vironment, and that is capable of autonomous 
actions in this environment in order to achieve 
its delegated objectives [5]. ABM is therefore 
inherently predestined to model and analyse 
decision makers as individuals and their inter- 
action and hence was chosen for this task in 
CONFIDENCE. 
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Figure 4 displays the concept of a software 
agent in more detail. The agent to the right in- 
teracts with the environment and also other 
agents on the left. It perceives observations on 
the environment by its sensors. The same way 
itis capable to perform actions on the environ- 
ment with its actuators. Agents are capable of 
evaluating their environment and of decision 
making by means of simple rules up to artificial 


intelligence. 
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Figure 4: Structure of an autonomous (software) agent in- 
teracting with its environment. 


It is obviously the decision making, that defines 
an agent in the end. In the context of modelling 
nuclear emergency management, an agent 
has to select the best countermeasure strategy 
from a set of strategies on the basis of param- 
eters. In CONFIDENCE we interviewed emer- 
gency managers from different countries by 
questionnaires, on how they make their deci- 
sion, what rules they follow, etc. From this in- 
formation a set of agent types was modelled 
and implemented as well as methods for nego- 
tiation between agents to find to a common 
consensus. An intelligent strategy evaluation 
system based on the agent-based negotiation 
simulation has been introduced in order to sim- 
ulate the decision making process of stake- 
holders on computationally tractable assump- 
tions. In the framework of the system, agents 
can score the recommended strategies before 
negotiation and negotiate them by using differ- 
ent negotiation skills [6,7]. Moreover, one indi- 
cator was introduced to describe how much the 
agents can compromise in each negotiation 
[7]. This parameter may also reflect the degree 
of selfishness behaviour of agents. Figure 5 
displays a chart of the demonstrator, where 
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agents negotiate on the ranking value of a 
strategy. 


Negotiation for one strategy 
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Figure 5: A group of 12 agents agreeing on a ranking 
value for a one strategy within 14 iterations. 


Combining methods 


Decision making methods have their own spe- 
cific upsides and downsides, e.g. an upside of 
MCDA is that it is simple and increases trans- 
parency, while as downside the alternatives to 
rank have to be known beforehand. Carefully 
combining decision making methods can 
greatly improve the overall usability as well as 
the quality of results. In the following two com- 
binations of methods are presented that can 
benefit from each other. 


Evolutionary Algorithms, ABM, and 
MCDA 


Evolutionary Algorithms (EA) are constructive 
optimisation methods, i.e. they search for an 
optimal, or in praxis close to optimal, solution 
within given parameter ranges, especially if no 
knowledge is available how to construct the 
optimal solution [8]. Regarding emergency 
management this means existing counter- 
measure strategies are modified in a way that 
they fit to the current emergency situation. 


The basic idea of EA is to encode a solution 
respectively countermeasure strategy as a ge- 


nome of an individual and to have a large pop- 
ulation of individuals evolve in the desired way. 
For this, each individual is evaluated in respect 
to its survivability, called fitness. The least fit 
individuals are removed and the population is 
filled up with new child individuals derived from 
the ones that survived by combining and mod- 
ifying the genes of some parents in a process 
called crossover and mutation. The process is 
repeated until some criterion to stop is 
reached. While the method is straightforward 
the challenge is obviously the encoding of a 
strategy as genome and the evaluation of the 
fitness of the individuals. 


In a complex decision support system based 
on ABM, e.g. for simulation of power supply 
management of a city, the EA can be intro- 
duced as a “super-agent” that constructs strat- 
egies and interacts with the agents for fitness 
evaluation. In addition to the negotiations with 
other agents, the agents as autonomous indi- 
viduals have their personal preferences on the 
suggested strategies. Therefore, it seems evi- 
dent to model such individual behaviour by 
MCDA. That way the advantages of the three 
methods can be combined into a more efficient 
decision support system. Currently we follow 
this approach in the framework of the HGF 
portfolio security for evaluating power distribu- 
tion management in future cities. 


Case Based Reasoning and MCDA 


As mentioned above, one integral and chal- 
lenging part of MCDA is to define alternative 
actions to be evaluated. An EA is one oppor- 
tunity to construct strategies to be further ana- 
lysed according to different criteria. Another 
approach is prepare a database with different 
scenarios and strategies beforehand to be 
used in a decision situation. This approach 
was particularly pursued in CONFIDENCE 
where Case Based Reasoning [9], a problem- 
solving paradigm, was applied to select strate- 
gies to be negotiated in the framework of ABM 
of the decision making process [6]. 
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Hence, for defining decision alternatives, CBR 
can be applied to limit the set of possible op- 
tions. Originating from cognitive science, CBR 
utilizes specific knowledge of previously expe- 
rienced problem situations to solve a new 
problem. The main assumption is that similar 
problems have similar solutions. Hence, a pos- 
sible set of decision alternatives is based on 
experiences of similar and solved problems, 
which in our context, correspond to simulation 
results. This approach and especially the de- 
termination of suitable criteria have been fur- 
ther elaborated in the framework of nuclear 
emergencies, where MCDA is applied to as- 
sess several possible disaster management 
strategies [10]. The objectives have been to 
rank different strategies in a transparent man- 
ner, to provide a broad discussion basis, and 
to preserve flexibility to account for the varia- 
bility of disasters and users' preferences. In 
particular, the contributions of the different cri- 
teria to the overall assessment are revealed. 


The basic idea for assessing strategies is to in- 
tegrate commonly discussed approaches that 
refer to performance measures and investigat- 
ing robustness. Furthermore, CBR related val- 
ues that reflect the trustworthiness of the solu- 
tions proposed are respected as well. The 


multi-criteria assessment considers current 
conditions, possible future developments, uti- 
lizes simulations of strategies to account for 
current constraints and uncertainties with re- 
gard to time, for example, and facilitates users' 
trust and understanding in the mechanism of 
the decision support method by integrating 
confidence values. These different perspec- 
tives are summarized (Figure 6) where the 
overall objective is to protect public and envi- 
ronment being decomposed into the criteria 
'effectiveness', ‘resources’, 'robustness' and 
'confidence'. Specifically oriented towards nu- 
clear emergencies, the effectiveness is meas- 
ured according to (i) the factor of dose reduc- 
tion (ii) the amount of waste, and (iii) public 
acceptance, taking into account non-radiologi- 
cal quantities as well. The criterion 'resources' 
states through which means the objectives are 
achieved and hence which resources and to 
what extent they are utilised. These results are 
gained by simulating the strategies considered 
and analysing them according to duration and 
uncertainties in respect of potential delays dur- 
ing the implementation possibly causing a re- 
duced resource utilization. Here, a strategy 
mo-del that is based on Petri nets is used al- 
lowing to capture combinations of measures, 
their order of implementation, the objects 
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Figure 6: Hierarchy of criteria for strategy assessment [10, Figure 6.6] 
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measures are targeted at, and resources [11]. 
Simulation results of JRodos [12] particularly 
provide costs of strategies. The criterion 'ro- 
bustness' considers uncertainties with regard 
to the extent of a disaster, changing environ- 
mental conditions, or insufficient information. 
For judging how robust a strategy is, a sce- 
nario-based approach [13] is pursued, investi- 
gating different scenarios and determining cor- 
responding effectiveness values under these 
varying conditions. The criterion ‘confidence’ is 
related to CBR and can be made measurable 
by similarity values and deviations in the differ- 
ent solutions. 


The multi-criteria assessment helps to struc- 
ture the decision problem, reduces its com- 
plexity, and promotes discussions of the stake- 
holders involved by, for instance, visual 
support (Figure 7). In particular, different views 
on strategy assessment are integrated ad- 
dressing various preferences that need to be 
respected in the final decision. Furthermore, 
different strategies can be discussed and ana- 
lysed according to their sensitivity in respect of 
weights (Figure 8) or criteria values. The strat- 
egy assessment particularly respects different 
temporal dimensions and hence current condi- 
tions as well as future uncertainties taking into 
account characteristics of the underlying deci- 
sion support method. 
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Figure 7: Assessment of strategies illustrated as stacked- 
bar chart depicting the contributions of criteria values to 
the overall assessment [10, Figure 6.8] 


» ua N 


Assessment 


| Strategy A — Strategy B 


Figure 8: Stability of result according to changes in the 
weight of factor of dose reduction’ [10, Figure 6.10] 


Summary and Future Work 


Multi Criteria Decision Analysis has proven to 
be a valuable tool for decision making in vari- 
ous interactive workshops. MCDA is relatively 
easy to apply once the criteria and preferences 
are agreed upon. It is fast in computation and 
therefore very interactive in its handling. The 
setup forces discussions among decision mak- 
ers and leads to clarity as well as transparency 
in the decision process. Considering uncertain- 
ties in the input improves the interpretation of 
the results and allows to determine more ro- 
bust solutions compared to the hitherto exist- 
ing approach. As a drawback applying MCDA 
requires knowledge on the method and the 
consequences of choosing specific criteria or 
normalisation methods. The setup takes some 
time as discussions are required among the 
decision makers to reach a consensus. There- 
fore, it is best suited for preparation or late 
phase emergency management where time is 
less pressing. 


Combining MCDA with other decision support- 
ing methods like Agent Based Modelling, Evo- 
lutionary Algorithms, or Case Base Reasoning 
can improve the performance of a decision 
support system, especially for complex prob- 
lems like the management of power distribu- 
tion in urban areas. The growing complexity of 
dependencies between critical infrastructures 
and the ongoing urban transformation towards 
smart cities, challenge crisis management. In 
particular, there is a lack of knowledge on pos- 
sible disruption scenarios, the range and se- 
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verity of cascading effects as well as appropri- 
ate management strategies. In the framework 
of the HGF portfolio security, we are working 
on robust and comprehensible solutions for cri- 
sis management, specifically for maintaining 
security of supply and protecting critical infra- 
structures in complex crisis situations. This es- 
pecially requires an understanding of emer- 
gences resulting from numerous interacting 
system components in an urban area. We 
have developed an agent based optimisation 
framework that will be further enhanced by, for 
example, multi criteria analysis capabilities for 
agents. Besides global strategies and objec- 
tives, individual agents aim at self-preservation 
and demand-driven supply of services, chal- 
lenging the assessment of potential strategies 
and opening up various research possibilities 
in the context of MCDA. 
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Introduction 


In 2019 the Hydrogen Group continued to de- 
velop the in-house specialized CFD codes, 
GASFLOW and COMSD, and conducted sev- 
eral experimental programs to deepen the un- 
derstanding of the mixing behavior and the 
transient combustion phenomena, flame accel- 
eration and deflagration-detonation-transition, 
of accidentally released hydrogen. Members of 
the group are actively transferring their insights 
and expertise into the respective IEA Hydro- 
gen Task 37, the standards committees of ISO 
(TC 197), CEN/CENELEC (TC 6) and in the 
European Hydrogen Safety Panel. Among 
other third party funded projects the Hydrogen 
Group has been conducting ambitious experi- 
mental programs for the European Fuel Cell 
and Hydrogen Joint Undertaking (FCH JU) 
funded projects PRESLHY - Pre-normative re- 
search for the safe use of liquid hydrogen - and 
HyTunnel-CS - Pre-normative research for 
safety of hydrogen driven vehicles and 
transport through tunnels and similar confined 
spaces. The report for 2019 will focus on some 
highlights provided by the Hydrogen Group in 
the field of hydrogen safety in tunnels and con- 
fined spaces. 


Hydrogen Fire Suppression in Traffic 
Tunnels 


The Hydrogen Group supported by the Pro- 
Science team has performed more than 600 
tests to investigate the effect of mechanical 
ventilation on hydrogen jet structure and dis- 
persion. The experiment was integrated in the 


huge test vessel V220, see Figure 1a,b. A pow- 
erful fan was installed in one of the two main 
ports of the vessel. Between the ports a 6 m 
long test domain was arranged to mimic a 
scaled down cut-out of a tunnel with an acci- 
dental release of hydrogen. Thirty different 
configurations were examined by changing the 
H2 mass flow rate from 1 to 5 g/s and the air 
flow velocity in the range 1.5- 5 m/s with co- 
aligned, opposite and tranversal direction of 
the ventilation flow in relation to the hydrogen 
jet direction. 


Camara/Light 


(reflected) screen floor 


Figure 1: a Experimental set-up of tunnel ventilation ex- 
periments in HYKA V220 test vessel; 
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Figure 1: b View from downstream on the ventilated tun- 
nel segment with central fan 


The aim of the experimental campaign was to 
determine the hazard distances as function of 
hydrogen mass flow rate, ventilation air flow 
velocity and relative release direction. It was 
found that for all examined configurations, the 
mechanical ventilation led to a considerable 
size reduction of the pre-mixed cloud defined 
by the flammable mixture. 


MOLF-ISJUNOI 


1gH,/s; 1mm. 


nozzle; no flow 
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1g H,/s; 1 mm nozzle; counter-flow ^5 m/s 


(b) 


Relevant tests have been selected for upload- 
ing of the corresponding result files on the Hy- 
Tunnel-CS open data repository. These data 
will be further used for code validations. 


Hydrogen Fire Suppression in Traffic 
Tunnels 


Fires are representing typical hazards in tun- 
nels and some devastating fire accidents 
where initiating the development of the Euro- 
pean tunnel safety directive [https://eur-lex.eu- 
tent/EN/TXT/PDF/?uri-CELEX:32004L00548f 
rom=en]. As hydrogen fires show some quite 
distinct characteristics and different behavior 
compared to conventional fires and the current 
version of the European tunnel safety directive 
does not address hydrogen as an alternative 
fuel, the applicability of conventional mitigation 
technology has to be evaluated. The suppres- 
sion effects of installed fire safety systems and 
the special accidental conditions have been 


co-flow 


(c) 
Figure 2: Measured concentrations in horizontal plane at height of release nozzle; (a) without ventilation, (b) with 
counterflow ventilation (c) with co-flow ventilation 
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studied experimentally and theoretically by nu- 
merical simulations. 


The CFD code Fire Dynamics Simulator (FDS) 
[https://www.nist.gov/services-resources/soft- 
ware/fds-and-smokeview], developed by the 
US NIST and widely used in the fire safety 
community, has been adopted to analyze the 
hydrogen fire generated from an accidental re- 
lease of hydrogen. The simulations account for 
the effects of water sprays and/or in oxygen 
starving conditions. The main variations con- 
cern the leaking mass flow rates of the hydro- 
gen source, the mass flow rates of water injec- 
tion, and /or different droplet sizes. 


Figure 3 shows the simplified 3D tunnel section 
which is modeled geometrically together with 
three cars. One of the cars releases its high 
pressure hydrogen inventory via a nozzle on 
the lower rear side, as a primary fire is as- 
sumed to open the temperature activated pres- 
sure release device (TPRD). The release noz- 
zle pointing backwards from the concerned 
car, is marked with red color in Figure 3. 


(b) Transverse view 


Figure 3: Numerical mesh of the simulated road tunnel 
segment with 3 cars, 

one of them (on the right hand side) releasing hydrogen 
through a nozzle (red dot) 


The gas temperature and gas compositions 
have been computed with and without water 
spray and with a variation of oxygen depletion 


in the concerned tunnel segment, i.e., less 
than 21 vol. 96 in normal air. The simulation 
cases are put together into the following table. 
The transient thermal state of the gas in the 
control volumes, e.g., temperature and steam 
fraction, have been calculated. The results are 
used to derive recommendations whether the 
environment in the tunnel is suitable for fire 
control, evacuation and rescue operations un- 
der the heat release rates of the hydrogen fire 
and for the selected water spray and oxygen 
depletion conditions. 


The following preliminary results are obtained 
by the FDS simulations. The temperatures of 
as shown in Figure 4 indicate a cooling effect 
of the water spray on the fire, although the tem- 
perature decrease is not that much due to the 
relatively low mass flow rate of water in that 
simulation case. It is interesting that for some 
cases the gas temperature with spray is higher 
than without spray (compare black and blue 
line in Figure 4). The reason for this surprising 
behavior is, that the hydrogen combustion is 
intensified by the turbulence and mixing effect 
induced by the spray droplet momentum. 
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Figure 4: Gas temperature changes in upstream caused 
by water spray 


The humidity (steam fraction) is also computed 
as it is a classical indicator for the impact on 
human beings. Obviously, the humidity at tun- 
nel exit is increased due to the operation of wa- 
ter spray. It means that some droplets must 
evaporate and remove tangible heat of hydro- 
gen fire. However, although this might have an 
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Table 1.1: Calculation cases of hydrogen fires suppressed by water sprays or oxygen depletions 


Water spray Oxygen depletion 
Small mass flow rate | Large mass flow rate Slight Medium | Serious 
of water of water starving | starving | starving 
Small Large Small Large of O2 of O2 of O2 
droplet droplet droplet droplet 
Small 
mass 1 3 5 7 a [e e 
flow rate 
of H2 
Large 
mass 2 4 6 8 b d f 
flow rate 
of H2 


impact on first responders strategies, there is 
practically no effect of the increased humidity 
on the reaction kinetics. Another observation 
is, that a higher water spray mass flux does not 
further improve the situation. 


So, it may be concluded that up to a certain 
limit a larger mass flow rate of spray brings bet- 
ter cooling effect and produces a higher steam 
fraction in the gas mixture. In certain conditions 
the spray may increase turbulence levels and 
promote faster combustion. The simulations 
with high degree of oxygen depletion show that 
even under extreme conditions hydrogen is al- 
ways burnt almost completely and no flamma- 
ble mixture is leaving the simulation domain. 
So a transfer of unburnt hydrogen into another 
segment or into the ventilation ducts is very un- 
likely. 


Consequence analysis of hydrogen 
explosions in tunnels 


In the more unlikely case of late ignition of re- 
leased and pre-mixed hydrogen in the tunnel, 
fire phenomena become less relevant and mit- 
igation of the potential explosion effects be- 
comes more important. Without proper mitiga- 
tion, the blast waves generated in a local 
hydrogen combustion event might travel a 
much longer distance in the tube like tunnel 
structures than in a free environment. To pro- 
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vide a reference scenario a typical tunnel acci- 
dent with late ignition was analysed with the in- 
house code GASFLOW-MPI [1]. This refer- 
ence scenario consisted of a tunnel segment 
with two parallel lanes in the same direction 
(two separate tubes), with a very mild inclina- 
tion, but without ventilation or spray activation, 
i.e. without accounting for any active mitigation 
technology. 


Figure 5: Tunnel segment modeled for late ignition / 
explosion scenario 


It was assumed that the inventory of 5 kg 
stored at one hydrogen driven car under high 
pressure is released via the TPRD at the rear 
of the car, similar as for the fire simulations 
above. For the release phase GASFLOW-MPI 
determines the time dependent mass of hydro- 
gen included in the flammable cloud and the 
potential for a deflagration-to-detonation tran- 
sition (DDT) in this premixed system (see Fig- 
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ure 6 top). This information was used to deter- 
mine the most conservative location and timing 
for an ignition. So ignition was assumed for 
maximum flammable mass at about 12.5 s af- 
ter start of the release in a central position of 
the premixed cloud with a hydrogen concentra- 
tion close to stoichiometry (about 3096). The ig- 
nition location close to the central ceiling may- 
be easily motivated, as in tunnels with such a 
horse-shoe shaped cross-section the electric 
supply lines for light and ventilation are instal- 
led typically there. However, an accidental ig- 
nition by a spark is considered a weak ignition 
Source, as it won't initiate a detonation directly. 


After the weak ignition GASFLOW-MPI deter- 
mines the transient pressure loads generated 
in the deflagration. Although a DDT may not be 
ruled out, because of the high degree of con- 
finement and obstruction, it is difficult to deter- 
mine the exact location and further develop- 
ments of a detonation. However, the loads 
generated by the deflagration, with local max- 
ima in the order of 3 bar (see Figure 7), are 
strong enough to damage all involved vehicles 
seriously and eventually to destroy the tunnel 
structure. 
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Figure 6: a): time evolution of flammable mass and 
relative diameter of DDT capable cloud; 
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b) Y-z-cut of hydrogen concentration and the location 
of ignition 


Figure 6: b) vertical cut of flammable cloud at critical 
ignition time (12.5 s) and critical ignition location (red 
dot) 


Figure 7: Blast wave caused by hydrogen deflagration 
reaches the car roofs 


It has to be stressed, that this scenario is as- 
sumed to have very low probability. It is mainly 
used as a reference for evaluating the effec- 
tiveness of mitigation measures against explo- 
sion scenarios. The next chapter describes 
such an evaluation for water sprays, which are 
typically installed in larger tunnels. 


Study of attenuation effect of water 
droplets on blast waves 


To gain a better understanding, first the differ- 
ent attenuation mechanisms, like momentum 
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absorption, thermal absorption, reflection and 
droplet breakup, were systematically assessed 
with simplified analytical models. Influencing 
parameters, such as droplet size and liquid 
phase concentration were analysed. A litera- 
ture survey provided an overview of similar re- 
search and allowed to identify high quality val- 
idation data generated in well instrumented 
experiments. 


A suitable droplet model, addressing the main 
phenomena droplet drag forces and droplet 
breakup, has been developed and imple- 
mented in the COMSD code. The correlations 
and parameters needed for the drag coefficient 
and for the breakup time and secondary drop- 
let size have been tuned to reproduce the re- 
sults published in [A. Chauvin et al. "Investiga- 
tion of the attenuation of a shock wave passing 
through a water spray”. In: International Sym- 
posium on Military Aspects of Blast and Shock 
(MABS) 21 (2010]. Figure 8 depicts the chan- 
nel used in these experiments with the high 
pressure driver region, a low pressure blast 
wave travel zone and the droplet cloud zone. 


The results achieved with the COMSD imple- 
mented model are compared to the experi- 
mental results in Figure 9. The simulations re- 
produce the results qualitatively and quanti- 
tatively reasonably well. In particular the early 
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critical phase of droplet atomization is captured 
very well. 


Figure 8: Experimental set-up for the Chauvin experi- 
ments 


Subsequently, using the developed model, 
simulations have been conducted to plan mid- 
dle sized water spray / shock wave interaction 
experiments in the hydrogen test facility V220 
of the hydrogen test center HYKA. Hydrogen is 
detonated in a combustion unit, what induces 
a shockwave. The pressure propagation in the 
absence and presence of a droplet cloud has 
been computed. The results will be used to 
tune the experimental layout. 


Finally, it is planned to use the COM3D code 
with these new extensions to assess the miti- 
gating effects of conventional water sprays and 
modern mist system for hydrogen explosion 
scenarios in traffic tunnels, as described 
above. 


— with Droplets (computations) 
=— no Droplets (computations) 
=-=- with Droplets (experiments) 
==- no Droplets (experiments) 
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Figure 9: Comparison of COMSD results and measured pressure at point 3190 without (red lines) and with droplet 


cloud (blue lines) 
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Study of a liquid air energy storage system 


Mathieu Champouillon, Elisabeth Schröder, Dietmar Kuhn 


With the emergence of renewable energies, 
the grid becomes more and more complicated 
to manage. The intermittency of these new 
means of production increases the need for 
energy storage solutions. 


Against this background the idea of an innova- 
tive energy storage system arises which uses 
electrical energy to liquefied air, then store the 
liquid phase in a tank in order to compensate 
the electricity surplus of the grid, and expand it 
through a turbine when electricity is needed on 
the network. 


The principle of function is given in the follow- 
ing diagram, Fig. 1. 


El Heat 


Heat 
: Liquid 
Ar zn g "mu EA» 
Storage Tank 
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we 
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Figure 1: Operating principle of an energy storage system 
based on air liquefaction. 


In the liquefaction part, we need electrical en- 
ergy to run the compressors and the pumps. 
During this process, thermal energy could be 
extracted from the compression (losses), and 
from the cooling of air to its boiling tempera- 
ture. For the expansion part, heat is needed to 
evaporate the liquid air and to improve the ex- 
pansion in the turbine. 


For this study, only the liquefaction part will be 
investigated. The feasibility and efficiency con- 
sideration of the evaporation system will be ex- 
amined in a future work. 


The liquefaction of air goes back to Carl von 
Linde and W. Hampson who both registered a 
patent for gas separation in 1895. It is the eas- 
iest cycle and contains one or two heat ex- 
changers for cooling the compressed air by the 
recirculating cold gaseous phase below the in- 
version temperature. 


Later, many other circuits were developed to 
improve energy efficiency through additional 
cooling and compression stages upstream of 
the expansion valve [1-9]. The Claude Cycle 
for instance uses three heat exchangers in a 
row. After the first heat exchanger, the com- 
pressed air flow is devided in two flows. One 
air stream passes through the liquefaction line, 
the second stream is expanded in a turbine 
and after mixing the turbine exhaust air with the 
recirculating cold air of the liquefaction line it is 
used as cooling media in heat exchanger two. 
Some of the compression energy can thus be 
recovered but additional plant components are 
necessary. Another variant of the Claude Cy- 
cle is the Kapitza Cycle, whereby the third heat 
exchanger is removed. Although both cycles 
are more efficient than the Linde-Hampson Cy- 
cle this study deals with the latter for simplicity. 
In order to find out if the air liquefaction process 
can be coupled to power generation pro- 
cesses, for instance ORC-Cycles, the influ- 
ence of parameter variations on a Linde- 
Hampson process is studied. 


In Fig. 2 the Linde-Hampson Cycle is illus- 
trated and the process steps are numbered. 
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Figure 2: Linde-Hampson cycle 


At ambient pressure, fresh air (1) is mixed with 
the recirculating cold air (9) to the inlet air (2) 
of the compressor (C). The mass flux at posi- 
tion (2) is m. After the compressor, the air pres- 
sure is enhanced to 200 bar (3). The com- 
pressed air is then cooled in an externally 
chilled heat exchanger to state (4). A second 
heat exchanger, operated by internal heat ex- 
change, reduces air temperature close to air 
condensation temperature (5). After isenthal- 
pic expansion in the throttle valve air pressure 
is reduced to ambient pressure (1.015 bar) 
thus causes the temperature to be lowered to 
the condensation temperature of -194.3?C (6). 
In the separator, gaseous (8) and liquid phase 
(7) are separated and the gaseous air (8) is fur- 
ther used in heat exchanger 2. 


The air states can be calculated as follows 
based on mass flux m at point 2, whereby x is 
the mass fraction of gaseous air: 


Ambient condition: T1 = 293.15K, pi = 
1.015 bar, fixed val- 


ues 


T; = xT + (1 — x)T,, 
p2 = pı, whereas x 
corresponds to 
mass fraction of re- 
circulating air 


Ideal mixing of fresh 
and recirculated air 


Isenthalpic expan- hs = he , ps = ps 
sion at valve: 
Boiling point atam- | Te = -194°C, pe = 
bient pressure: 1.015 bar 
Gaseous phase at T7=Te, pz- pe, 
boiling point m; —-(1—x):m 
Liquid phase: Ts= Te, pe- pe 
, Mg = mMm'x 


Adiabatic cooling in 
HE 2 


Tə = Ts + (T4 -— T5) / 
X, po = pe 


Isothermal compres- 
sion: 


Ts = T1, ps = 200 
bar, fixed values 


Cooling in HE 1: ha = (1 — Neoo1) * ha 
with efficient coeffi- 
cient Noo) T4 can be 
calculated from h4 
by using REFPROP 
database for in- 


stance, p4 = ps 
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What can be seen from the formulas above, x 
and Nco0ı cannot be varied independently. 
Some physical limitaions have to be consid- 
ered. In both heat exchangers, the cold fluid 
cannot be hotter than the hot one. That leads 
to following restrictions in HE 2: T,» 


To and T; > Tg. This could not be respected 
when x Or noo are too low. In Fig. 3 tempera- 
tures T4 and Ts are shown as function of nego) 
for given x=0.5. 


Evolution of T4 end TV in fonction of the cooling coefficient for x 0.5 and T3 = 293.15 K 


Figure 3: Temperatures T4 and Ts depending on 7.001 for 
given x = 0.5 


As shown in Fig. 3 for x = 0.5, Ncocı has to be 
higher than 0.74 which means that the effi- 
ciency of HE 1 is restricted in order to provide 
a minimum cooling rate. Depending on the 
amount of liquid air (1-x), the first stage cooling 
in HE 1 is predefined. In addition following cri- 
teria for HE 1 has to be fulfilled as well: T, > 
Tio and T; > T,,. Therefore, a specific external 
cooling process must be adapted to the spe- 
cific temperature level. 


In order to understand the behaviour of the 
system, the influence of the liquid rate (1-x) 
and the cooling coefficient Ncocı on the temper- 
atures and the enthalpies of each point are an- 
alysed. The interesting values are: 
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Figure 4: Multi-stage compression with heat extraction 


e Temperature at point 4 (T4) 
e Temperature at point 9 (To) 


e  Enthalpy extracted between 3 and 4 
(Aha4) in HE 1 


e Maximal temperature of the compressor, 
like if there were not any cooling (Tcomp) 


e  Enthalpy that could be extracted by the 
compressor cooling system 


The respective parameters are varied in order 
to examine how much heat can be extracted 
by combining air liquefaction with external pro- 
cesses. Possible combinations are identified. 
First, compression heat is used for warming up 
the working fluid of an external ORC-process 
of a geothermal energy plant. For this purpose 
the air compression must be operated in sev- 
eral stages in order not to overheat the com- 
pressor. Second, the influence of liquid air ratio 
(1-x) on the efficiency of HE 1 (Ncocı) and third, 
the extraction of gaseous air prior to mixing are 
investigated. 


Influence of parameters 


In order to study the influence of parameter 
variation like x or Ncooı the isothermal compres- 
sion is supplemented by a multi-stage adia- 
batic compression for both enhancing the total 
energy extraction from the circuit and avoiding 
overheating of the compressor during adia- 
batic operation, as shown in Fig. 4. This 
method allows to increase the released energy 


yield and to use it technically in another pro- 
cess. 


For each stage, we consider an adiabatic com- 
pression and then an isobaric heat exchange. 
The heat exchange can be effected by another 
system (ORC p.e.). The compression ratio is 
the same for each compressor: 


za — constant at each stage. (1) 
It is important to notice that the temperature af- 
ter each heat exchanger (including Ts) de- 
pends on the system used to cool down the air, 
because its minimal value is the temperature 
of the inlet of the cooling fluid. After the last 
stage, it is assumed, that T3 remains at 20°C 
like in the isothermal case. The total heat which 
can be used in external processes is defined 
as heat extracted in the compression stages 
and heat released in HE 1. As described below 
the total heat release depends on x and Neoo 
and is shown in Fig. 5. 


Influence of Ncoo on heat extraction 

in HE 1 

In this examination the liquid air ratio is kept 
constant at (1-x) = 0.1 whereas Ncooı is varied 
between 0 and 0.7. The efficiency of HE 1 is 
determined by the value of Ncooı whereby Ncooi = 
0 indicates that no heat is extracted in HE 1. 
Both a one-stage and two-stage compression 
is regarded and the total heat, provided by ad- 
iabatic compression and by HE 1 is shown in 
Fig. 5 
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Figure 5: Evolution of the total recovered enthalpy for 
different values of Noo 
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Here we can see that when Nooo is increasing, 
the total enthalpy recovered by cooling sys- 
tems Ahiota is decreasing. While the enthalpy 
between 3 and 4 that is extracted in the HE 1 
is higher, the compression enthalpy is even 
lower because T2 is decreasing a lot, until 111 
K for Nooo = 0.7. 


If Ncooı were greater than 0.7, the temperature 
Tg would be lower than the boiling temperature, 
which is impossible because we need gaseous 
air in this part. 


When the number of compression stages in- 
creases, the total enthalpy recovered by cool- 
ing systems is lower (-36% for 1 to 2 stages), 
because the maximum temperature of the 
compression is lower. But the influence of Neooı 
is also reduced: for a one-stage compression, 
the difference of total enthalpy recovered is 
230 kJ/kg for 0 < Nooo < 0.7, but only 80 kJ/kg 
for a two-stage compression. 


The choice of Ncocı depends on the type of cool- 
ing system. If the user prefers to give a lot of 
energy through HE1, he will have better to in- 
crease Ncool, but if he wants to recover as much 
energy as possible, he will have to decrease 
Ncocı to take advantage of the compression. It 
has to be mentioned here, that a one-stage ad- 
iabatic compression leads to unreasonable 
high temperatures of 773°C which cannot be 
tolerated because which would lead to the de- 
struction of the compressor. Taking a temper- 


60 


ature limit of 500°C for the compressor into ac- 
count, a one-stage compression cannot be re- 
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Figure 6: Total enthalpy recovered as function of x and 
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As a consequence of the increasing yield of liq- 
uid air, the recirculated cold air (indicated as x) 
decreases, whereby the mixture temperature 
T2 is increased. Therefore, more energy is re- 
leased, in the compression heat exchangers 
and in HE 1, so that the total enthalpy 
(delta H total) increases, as shown in Figure 
6. One more time, it is important to notice that 
some values are physically impossible: Ncooı < 
0.4 for x = 0.6 and Ncooı > 0.7 for x = 0.9. 


The goal of this storage system is to liquefy air 
in order to store energy. When x decreases, 
more liquid air can be produced. The previous 
results are therefore good news since the sys- 
tem allows to create more liquid air, while re- 
covering more energy in the cooling systems. 
But the main problem of decreasing x is the 
significant limitation of the operating ranges for 
many parameters. For example, Ncocı must be 
between 0.4 and 0.8 for x=0.6. 


Influence of temperature T3 (entrance 
of HE 1) 


This temperature corresponds to the minimum 
temperature of the compressor cooling system 
and the maximal temperature for the HE 1. 
This means that the value of Ts will be decisive 
concerning the choice of the system used to 
cool down the air in the compressor and the 
first heat exchanger. 
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Comparison of the influence of the compressor and the heat 
exchanger 


—®— Compression T3 = 20°C 
—#— Compression T3 = 60°C 
—®— Compression T3 = 120°C 
—9— HE1 T3= 20°C 
—9— HEL T3 = 60°C 


—9— HE1 T3 = 120°C 


Figure 7: Influence of T3 on recovered enthalpy from 
compressor cooling system and HE 1 


The temperature of 60°C corresponds to the 
minimum temperature of a district heating sys- 
tem, and the temperature of 120°C has been 
chosen arbitrarily to see the behaviour of the 
enthalpy at high temperature. It can be seen 
from Fig. 7 that the higher the temperature Ts, 
the higher the total enthalpy recovered in the 
compressor and in the HE1. The second effect 
of the increase of temperature is the tightening 
of the range of Ncooı values, for instance: with 
Ts = 120°C, ncooi must be higher than 0.2. From 
Fig. 7 it can further be seen the impact of the 
enthalpy recovered during the compression. 
For low values of Ncooı, the difference of the to- 
tal enthalpy recovered is mainly due to the 
compression, with increasing impact of HE 1 
for higher values of Neooı. 


Influence of cold air extraction between 
point 8 and 9 


At the point 9, we have got cold air with a mass 
flow rate of x that could be used for another 
application instead of being recycled in the cy- 
cle at the point 2. The idea is to extract a cer- 
tain amount of air at the point 9, with a mass 
flow rate of y. This modification has not a sig- 
nificant impact on the system because it modi- 
fies only the temperature 2, which tend to Tam. 
bient When y tend to x. 


Enthalpy recovered during the compression and in the HE1 (x = 0.8) 


(ki/kg) 


Total enthalpy 


Figure 8: Influence of air extraction ratio y prior to mixing 
on total recovered enthalpy 


The main result as can be seen from Fig. 8, is 
that the enthalpy extracted during the com- 
pression and in the HE1 is increasing with y 
and ncooi. This happens because the tempera- 
ture T2 is increasing when the remaining mass 
flow rate of the point 9 decreases, so the en- 
ergy released by the compression will be 
higher. 


Conclusion 


This study provides a more detailed insight into 
the Linde-Hampson cycle for the liquefaction of 
air as an energy storage medium. Many pa- 
rameters were varied to investigate their influ- 
ence on the overall process. As a result of this 
study, coupling possibilities to other power 
generation processes, e.g. ORC cycles, are 
easier to identify. 
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Blockage formation experiments in a water rod bundle 


Karsten Litfin, Jasmin Klupp, Andreas Class 


Introduction 


Solid particles flowing into (or within) the core 
of a nuclear reactor can block one or more sub- 
channels of a fuel assembly. As they do not af- 
fect substantially the outlet conditions (temper- 
ature, pressure), they can go unnoticed by the 
global instrumentation. However, they lead to 
local heat transfer degradation and the corre- 
sponding increase of the cladding tempera- 
ture. Under certain conditions, depending 
among others on the size, location and poros- 
ity of the blockage, this can result in local fuel 
cladding failure. Thus, a deep understanding of 
the formation and main characteristics of these 
blockages is essential for the development and 
safety assessment of innovative reactor con- 
cepts based on heavy liquid metal coolant, i.e. 
with particles, that have usually lower density 
than the coolant. 


Former intensive test programs on particle for- 
mation in sodium and water-based systems 
[1], [2] cannot or can only limited be transferred 
to heavy liquid metal (HLM) cooled reactors, as 
they are based on particles with higher density 
than the coolant [3]. For this reason, dedicated 
experiments with the proper density ratio are 
mandatory. 


For these tests, a 19-pin hexagonal fuel bundle 
mock-up in a transparent flow channel sub- 
jected to a water flow seeded with buoyant par- 
ticles of relevant density, size, shape and wet- 
ting properties has been set up at the 
Karlsruhe Liquid Metal Laboratory correspond- 
ing to the dimensions of the Lead-Bismuth 
(LBE) cooled MYRRHA research reactor pro- 
jected in Mol, Belgium. 


The nature of the blockage as well as the for- 
mation is characterized by optical measure- 
ments, supplemented by flow rate, pressure 


drop and temperature measurements. For the 
handling of the particles, both an injection (up- 
stream of the test section) and extraction 
(downstream) subsystem has been imple- 
mented into the existing facility. The density ra- 
tio liquid/particle was to fulfil the applicable 
scaling laws between the water tests and the 
HLM reactor. 


Experimental Setup 


The KALLA water loop was originally used for 
water rod bundle experiments in the framework 
of EU-Project DEMETRA. A complete descrip- 
tion of the experiments and the loop design can 
be found here [4]. 


The maximum flow rate of the loop is 130m3/h 
with a pressure head of 14.7 bar. The active 
loop inventory is 801, the maximum total inven- 
tory is about 8m?. The loop is equipped with a 
heat exchanger to keep a constant fluid tem- 
perature and a cleaning system to remove 
used particles. Due to the PMMA test section 
the temperature is limited to 50?C. The flow 
rate is measured by a commercial flow-meter. 


Rod bundle test section 


The geometrical details of the 19-pin rod bun- 
dle for the blockage formation and growth ex- 
periment is shown in Figure 1, the comparison 
between the experimental setup and the con- 
ditions in the MYRRHA reactor rod bundle in 
Table 1. A detailed figure and a picture of the 
test section including the particle injection is 
given in Figure 2. 
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Figure 1: Water rod bundle geometry 


Table 1: Geometric parameters used for the test section 
compared with MYRRHA reactor core geometry 


MYRTE test 


Parameter Symbol at KALLA MYRRHA Ratio 
Number of rods, - N 19 127 0.15 
Rod outer diameter D, mm 6.55 6.55 1 
Distance between rod centers (pitch) P, mm 8.378 8.378 1 
Pitch-to-diameter ratio P/D, - 1.279 1.279 1 
Flat-to-flat distance FF,mm 395 97.55 0.4049 
Wall distance W, mm 8.5139 8.5166 zl 
Wall distance-to-diameter ratio W/D, - 1.2998 1.3002 zl 
Wire diameter d, mm 1.80 1.80 1 
Wire pitch H, mm 262 262 1 
Wire pitch-to-diameter ratio H/D, - 40 40 1 
Total length L, mm 655 1400 0.467 
Total length-to-wire pitch ratio L/H, - 25 0.467 
a (complete bundle) Awr mm? 662.6 0.1821 
liameter (complete bundle) dhs mm 4.173 1.052 
ea (internal sub-channels) Asch, mm? 12.273 1 
Hydraulic diameter (internal dysch, mm 3.743 1 


sub-channels) 


(1:2)4 


oe 


Figure 2a: Detailed drawing of the foot part of the test 
section with the particle insertion (left) and the optical 
parts of the MYRTE test section (right). 
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Figure 2b: Picture of the mounted test section 
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Instrumentation 


The PMMA housing offers direct optical access 
to the blockages through the PMMA casing of 
the rod bundle. However, this access is re- 
stricted to the outer sub channels only as 
shown in Figure 3. 


Figure 3: Picture of a blockage in the wall sub channel 
through the PMMA casing 


The blockages of the inner sub channels are 
optically detected by an endoscope that is in- 
serted in one of the PMMA tubes of the rod 
bundle as depicted in Figure 4. As the endo- 
scope is equipped with it's own light source, 
further illumination of the investigated part of 
the rod bundle is not necessary. 


The use of this endoscopic system is a result 
of a comparison of different measurement 
techniques for the detection of the expected 
blockages by J. Biernath in his Master thesis 
[5]. It turned out that acoustic measurements 
are not suitable for this experiment as the com- 
plicated movement of the UDV (ultrasonic dop- 
pler velocimetry) sensor heads and the ex- 
pected interference signals prohibit a reliable 
detection of the blockage. From the optical 
measurements the endoscope was the easiest 
and most reliable option. 


"S 

Plexiglass 

-Blockage 
-— Endoscope 

- Water 


Figure 4: Scheme of an endoscope inserted into the 
rod-bundle mockup for blockage investigation (top) and 
a detected particle (bottom) 


Experiments 


For the blockage experiments, polypropylene 
particles with the density of 0.946g/cm? were 
used to simulate possible blockage formation 
of Lead oxide in LBE as the density ratio is very 
similar as shown in Table 2. 


Table 2: Density of blockage material for the experiment 


Material Observation — PP MYRRHA — Pp.MYRRHA Pp,MYRTE 
kg mr PLE g cm“ 
Stainless steel (DIN 1.4571) At 2 
Le le (PbO) At 2 
Ire > (magne hts , Fe30,) At 2 
Oxide fue :1 (MOX, 


> [20] 7756 0.7474 0.746 
7] 9530 0.9184 0.917 
[7] 5180 0.4992 0.498 
% porosity) At 20°C [3] 10565 1.0181 1.016 


In total 4 different experimental series have 
been performed with particle sizes from 1.5mm 
to 4.0mm diameter: 
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e particles that cannot pass the pin fixer fol- 
lowed by smaller particles 


e particles that pass the pin fixer but are 
likely to get stuck in the sub channel 


e high amount of irregular shaped particles 
that get stuck at the pin fixer and cause a 
noticeable pressure loss 


e one pre blocked sub channel and small 
particles that should not cause a blockage 
formation. 


The first experiment of the series is presented 
in this report where in 3 steps 180 particles with 
3.0mm diameter followed by 70 particles with 
2.0mm diameter were inserted into the flow up- 
stream the rod bundle test section. 


Insertion of 20 particles with 3.0mm 
diameter 


Figure 5 shows the test section after insertion 
of the first 20 particles (black spheres). Most 
particles are captured in the pin-fixer as the di- 
ameter is too big for the inner sub-channels. 
One particle passed the pin fixer and finally 
was captured in a wall sub channel. During the 
20 minutes of insertion the flow rate decreased 
by 396 and the pressure loss of the test section 
increased by 196. None of the particles was 
able to pass the test section. 


Figure 5a: 19 particles of 3.0mm diameter in the pin 
fixer 
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Figure 5b: One particle in a wall sub channel 


Insertion of 160 more particles with 
3.0mm diameter 


During the second run, 160 more particles 
(white spheres) of 3.0mm diameter are in- 
serted into the primary flow. Again, most parti- 
cles are captured in the pin-fixer and a small 
fraction of 15 particles go into the wall sub 
channels. In total 3 wall sub channels are af- 
fected and blockages with an axial extension 
of up to 9 particles are visible as shown in Fig- 
ure 6. Further decrease of flowrate and in- 
crease of pressure loss is depicted in Table 3. 


FA 


Figure 6a: Insertion of 160 more 3.0mm particles. Pin 
fixer 
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Figure 6b: Insertion of 160 more 3.0mm particles. Af- 
fected wall sub channels 


Insertion of 70 particles with 2.0mm 
diameter 


After the pin fixer and the wall sub channels 
have been blocked with a total of 180 particles 
with 3.0mm diameter, the blockage was 
probed with the insertion of particles with 
2.0 mm diameter. 


Most of the smaller particles are captured in 
the pin-fixer by the present 3.0mm particles. 4 
particles moved into a pre-blocked wall sub 
channel and 1 particle was able to enter an in- 
ner sub channel and got blocked there. This is 
shown in Figure 7. During the complete run, a 
further flow rate decrease and pressure loss in- 
crease was measureable and is shown in Ta- 
ble 3. After the run, the complete test section 
was dismounted, cleaned and investigated. No 
particle was able to pass the test section. The 


set up of the experiment as well as the com- 
plete 1st run was part of a Bachelor Thesis and 
is described more detailed elsewhere [6]. The 
other experimental runs are also described in 
detail elsewhere [7]. 


Figure 7: Insertion of 2.0mm particles in green into the 
pre-blocked test section. Pin fixer (top), wall sub chan- 
nels (middle) and inner sub channel (bottom) 
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Table 3: Flow rate and pressure loss during the first ex- 
periment 


Particles |20@ | 1609 |20@ |» 
3.0 3.0 2.0 

Flow rate | 3,596 | 19,5% | 2,596 25,596 

decrease 

Pressure 1,596 | 996 2,596 1396 

loss in- 

crease 


Concluding remarks 


The 4 performed experimental series show 
some predictable but also some new results 
and insights into the blockage formation in wire 
wrapped rod bundles. 


The first experiment, a blockage using 3.0mm 
particles that attach at the pin fixer and the wall 
sub channel followed by smaller 2.0mm parti- 
cles shows that even with 180 particles not all 
sub channels are blocked and still smaller par- 
ticles can pass by. However, they are very 
likely to attach to the existing big particles. The 
changes in pressure loss and flow rate are 
even for this large amount of particles very low. 


The second experiment was a blockage using 
2.0mm particles that theoretically should be 
able to pass the pin fixer and the test section. 
In contrast to the geometrical preconditions 
several blockages with almost half of the in- 
serted particles was observed. 


The third experiment was a blockage using a 
huge amount of particles that mainly blocked 
the pin fixer but also formed large axial ex- 
tended blockages in the wall sub channels. 
Here, a high influence on flow rate and pres- 
sure drop was detected. 


The fourth experiment was the use of small 
1.5mm particles that easily should pass the pin 
fixer and the test section but could attach to an 
existing blockage. Here, it was observed that 
on the one hand they form temporary block- 
ages that disappear and on the other hand, 
that no particle attaches to the existing block- 
age even as the particles were inserted into the 
same sub channel. 
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Finally we can conclude from the experiments: 


e Most of the measured blockages have a 
quite small axial extension. 
The probability of blockage agglomeration 
with a large axial extension is not very high 
as the particle size must be large enough 
not to be able to follow the enhanced 
crossflow upstream the blockage on the 
one hand but also must be small enough 
to meet the geometrical restraints at the 
pin fixer to enter the sub channel. 


e All particle sizes from 1.5mm form block- 
ages 


e Existing Blockages enhance crossflow 
and reduce agglomeration of small parti- 
cles 


e 1.5mm particles are able to follow the sub 
channel crossflow, the 2.0mm particles 
aren't and agglomerate on a excising 
blockage. 
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